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ABSTRACT 


This  report  addresses  the  prediction  and  analysis  of  solar  eclipse  circum¬ 
stances  of  interest  to  atmospheric,  ionospheric,  and  solar  studies.  In  the  first 
section,  specific  algorithms  for  use  in  pre-eclipse  planning  and  post-eclipse  analysis 
are  presented.  In  the  second  section,  appendices  present  1  he  calculation  of  solar 
and  lunar  ephemeridcs  of  requisite  accuracy  for  prediction  purposes,  the  calcula¬ 
tion  of  ephemeris  sidereal  time,  the  prediction  of  the  shadow  outline  on  and 
above  *he  earth,  the  calculation  of  local  eclipse  circumstances,  and  the  develop¬ 
ment  of  solar  coordinate  systems,  to  provide  the  background  for  and  the  founda¬ 
tion  of  the  preceding  algorithms. 
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1.0  INTRODUCTION 


This  report  addresses  the  prediction  and  analysis  of  solar  eclipse  circum¬ 
stances  of  interest  to  atmospheric,  ionospheric,  and  solar  studies.  In  the  first 
section,  specific  algorithms  for  use  in  pre-eclipse  planning  and  post -eclipse  analysis 
are  presented  in  a  form  suitable  for  manual  or  machine  calculation.  In  the  second 
section,  five  appendices  provide  the  background  for  and  the  foundation  of  these 
algorithms. 

Although  much  of  the  material  in  this  report  is  familiar  to  the  astronomical 
community,  its  dispersion  in  the  literature  or  its  rendition  in  a  variety  of  forms  do 
not  make  it  immediately  useful  to  the  non-astronomer.  This  report’s  first  section 
is,  therefore,  a  compilation  of  familiar  results;  the  second  section  is  tutorial  in 
nature. 

Appendix  A  presents  algorithms  for  calculating  the  geocentric  coordinates  of 
the  sun  and  moon.  In  the  solar  case,  the  algorithms  are  based  primarily  upon  the 
Newcomb  Tables  [  7 1  and  secondarily  upon  the  more  recently  published  Tables  of 
Jean  Mecus  [51.  In  the  lunar  case,  the  algorithms  are  based  primarily  upon  the 
“Improved  Lunar  Ephemeris’’  (3)  and  secondarily  upon  Meeus.  The  series 
developments  of  solar  and  lunar  longitude,  latitude,  and  distance  (parallax)  have 
been  truncated  in  view  of  the  relaxed  accuracy  requirements.  In  all  other  respects, 
however,  the  developments  of  the  solar  and  lunar  ephemerides  are  identical  to 
those  given  in  the  fundamental  references. 

Appendix  B  outlines  the  calculation  of  Ephemeris  Sidereal  Time. 

Appendix  C  develops  algorithms  necessary  to  predict  the  outline  and  motion 
of  the  lunar  shadow  on  an  earth  spheroid  of  arbitrary  radius  and  flattening.  It  is 
based  upon  the  exposition  of  Chauvenet  (1]  and  the  authoritative  summary  in 
The  Explanatory  Supplement.  .  ,[2| .  It  is  repeated  here  because  the  summary  of 
requisite  formula  in  [2) ,  with  a  different  ordering  from  that  of  (1),  lacks  both 
the  justification  and  seeming  inevitability  with  which  elements  related  to  the 
flattened  spheroid  and  the  shadow  motion  thereon  were  originally  introduced. 
Furthermore,  it  seems  appropriate  to  show  that  the  relationships  pertaining  to 
ionospheric  center  line  position  and  duration  follow  naturally  from  a  develop¬ 
ment  in  which  scale  factors  are  explicitly  employed. 

Appendix  D  develops  the  algorithms  required  to  predict  all  circumstances  of 
an  eclipse  at  a  locale  once  the  coordinates  of  the  locale  have  been  specified,  to 
adjust  these  circumstances  by  means  of  differential  correction  procedures  for 
modest  departures  from  the  locale  specified,  and  to  calculate  other  topoccntric 
parameters  of  interest. 

Finally,  Appendix  E  describes  various  coordinate  systems  applicable  to  solar 
astronomy  and  defines  their  interrelationships.  Particular  attention  is  directed  to 
'he  development  of  “pointing  instructions”  for  observation  of  phenomena  both 
on  and  above  the  solar  surface. 
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2.0  PRELIMINARY  REMARKS 


l-phcmeris  Time  (ITT.),  which  is  related  to  Universal  Time  <U  T.)  by  means 
of  t lie  expression 


t.T.  =  U.T.  +  AT  (I) 

will  be  employed  throughout  this  report,  except  where  indicated.  Similarly,  the 
ephcmeris  longitude  of  a  location  Xc,  which  is  related  to  the  geodetic  longitude  X 
(West  taken  as  positive)  by  means  of  the  expression 

\  =  \  +  1 .002738  AT.  (2) 

will  also  be  utilized,  except  where  indicated. 

In  either  instance,  conversion  irom  one  “system”  to  the  other  requires  the 
explicit  assignment  of  a  value  for  AT.  For  the  past  years,  values  of  AT  arc  given  in 
The  American  Eplu'mens  and  Sautical  Almanac,  page  vii,  luhle  <>)  Time- 
Difference  AT.  l  or  the  1970’s  Table  1  -  extracted  from  Tabic  <i7a,  Reduction 
from  Universal  to  Ephemeris  Time  of  1 21  -  provides  a  useful  extrapolation. 


TABLE  1 

AT  VALUES  FOR  1968-1980  INCLUSIVE 


1969.5 

+  42* 

19755 

+  51’ 

19705 

44’ 

1976  5 

53’ 

1971.5 

45‘ 

1977.5 

54’ 

1972.5 

47’ 

1978.5 

56’ 

1973.5 

4tj5 

1979.5 

57' 

1974  5 

501 

1980.5 

59’ 
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3.0  CALCULATION  OF  BESSELIAN  AND  AUXILIARY  ELEMENTS 


This  section  will  follow  the  development  presented  in  Appendix  C  and  is 
based  on  the  assumption  that  the  following  quantities  arc  available  at  the  time(s) 
of  interest : 


«0  $0 

R 

ac  >  4 


W 


r 


E.S.T. 


apparent  right  ascension  and  declination  of  sun 
lef  (A2!)| 

the  solar  radius  vector  in  A.U.  |cf  (A14)] 

apparent  right  ascension  and  decimation  of  moon 
[cf(A46)l 

the  horizontal  equatorial  parallax  of  the  moon 
(cf  ( A45 )  1 

Ephemcris  Sidereal  Time  |  cf  ( B3)  | . 


3.1  SHADOW  AXISPARAMETERS 


Calculate  the  shadow  axis  right  ascension  a  and  declination  d  and  the  solar- 
lunar  separation  g  (in  terms  of  R)  from 


where 


g  cos  d  cos  a 

=  cos  §0  cos  Oq  -  b  cos  6t  cos  a. 

(3) 

g  cos  d  sin  a 

=  cos  50  sin  Oq  -  b  cos  6,  sin  a. 

(4) 

g  sin  d 

=  sin  -  b  sinj, 

<5) 

b 

_  0.000  04  2  664 

R  sin  rr. 

(6) 

3.2  ''FUNDAMENTAL"  RECTANGULAR  COORDINATES  OF  MOON 


Calculate  the  rectangular  coordinates  x,  y.  z  of  the  moon  with  respect  to  the 
fundamental  plane(in  unitsof  the  earth’s  equatorial  radius  a  =  6378.160  km) 
from 


x  =  r. 

Icos  5,  sin  (a,  -  a)] 

(7) 

y  =  r* 

(  sin  6,  cos  d  -  cos  6,  sin  d  cos  (a,  -  a)) 

(8) 

z  -  rt 

I sin  6,  sin  d  +  cos  6,  cos  d  cos  (a,  -  a)] 

(9) 

5 


where 


rc  =  1  / sin 


(10) 


3.3  EPHEMERIS  HOUR  ANGLE  OF  SHADOW  AXIS 
Calculate  n,  the  cphcmeris  hour  angle  of  the  shadow  axis,  from 

#i  =  E.S.T.  -  a  (11) 

3.4  SHADOW  CONE  GENERATORS 

Calculate  /, ,  the  generator  of  the  penumbral  cone,  and  / j,  the  generator  of 
the  umbral  cone,  from  * 


sin/, 

=  0.0046  6401 8/gR 

(12) 

sin/j 

=  0.0046  40783/g  R 

(13) 

3.5  SHADOW  CONE  VERTEX  DISTANCES 

Calculate  c,,the  distance  of  the  penumbral  cone  vertex  above  the  fundamen¬ 
tal  plane,  and  c3 ,  the  distance  of  the  umbral  cone  vertex  above  the  fundamental 
plane,  from 


c,  =  z  +  0.2724  880cosec/,  (14) 

c2  =  z  -  0.2724  880  cosec /,.  (15) 

3.6  SHADOW  RADII  ON  FUNDAMENTAL  PLANE 


Calculate  C, ,  the  penumbral  radius  on  the  fundamental  plane,  and  the 


.  *  '  ,  -  J  * . .  - - *L-  C. J _ 

UJIIUIdl  IdUIUS  uil  iltc  ilUtUdiilWIiuit  puuiv,  iioiit 


*  tf ,  n  in  tne  discussion  leading  to  (E46),  we  wish  to  consider  the  radio  sun  whose  radius  is  3CR 
times  the  optical  radius  where  JCR“  (1  +  5CR),  then 

=  (0.0046  64018  +  0.0046  52367XR)/gR 
and 

sin/;,R  -  (0.0046  40783  +  0.0046  52367  JCR)/ g  R. 
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c,  tan  /, 


(16) 


C,  s 

Cj  ~  c2  tan  / 2  (17) 

where  C,  >  0;  C2  <0  for  totality,  i’2  >  0  for  annularity 

3.7  HOURLY  VARIATIONS 

Calculate  the  hourly  variationsd,  x,  y,p,  6,  and  S2  by  appropriate  numerical 
differentiation  of  the  Be&sclian  elements  sin  d,  cos  d,  x,  y,  p,  S|  and  t2 . 

3.8  AUXILIARY  ELEMENTS 


Cahulate  the  auxiliary  elements  P, ,  P2,  sin  d 
cos  (d,  -  d2)  from 

, ,  sin  d2 ,  sin  (d  ,  -  d2 )  and 

P,  =  (1  -e2  cos’d)* 

(18) 

p2  =  (1  -  e’sin2  d)* 

(19) 

sind,  =  sind /P, 

(20) 

cosd,  =  (1  -  e2)*  cos  d/P, 

C) 

sin  (d,  -d2)  =  e2  sin  d  cos  d/P,  P2 

(22) 

cos  (d,  —  d2 )  =  (1  -e2)W  /p,  p2 

(23) 

wnere  the  ellipticity  e  is  given  by  e  =  (0.0066  9454)* 

for  the  earth  spheroid 

Also,  calculate  the  additional  auxiliary  elements 

3. ,  a; .  b,  c.  and  c-,  from 

a,  =  _  j,  _p  xtan/,  cosd 

(24) 

a3  =  -  C2  -  p  x  tan  /2  cos  d 

(25' 

b  =  -  y  +p  x  sin  d 

(26) 

c,  =  x  +  pysmd+pS,  tan  /,  cos  d 

(27) 

c2  =  x  +  p  y  sin  d  +  p  £2  tan  /2  cos  d 

(28) 
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4.0  CALCULATION  OF  IONOSPHERIC  CENTRAL  LINE,  DURATION. 
AND  SHADOW  OUTLINE 

TTic  following  development.  except  when*  indicated. is  based  upon  Sections  (.'.9 
and  t  .  I  I  of  Appendix  C. 

4.1  IONOSPHERIC  CENTRAL  LINE  POINT 

I  rom  the  relation''. 

y,  =  y  Ip,  (29) 

f,  =  ITT-  -  x3  -  yj  )-‘  .  (30) 

where 

3<  =  I  +  0.15678503  x  10  J  h(  (31) 

in  which  h,  is  the  ionospheric*  center  line  height  above  sea  level  in  kilometers,*  cal¬ 
culate  a  point  on  the  central  line  from 

tan  (h1  = _ _  (32) 

J,  cos  d,  -  y,  sin  d, 

and 

f ,  sin  d,  <  y,  cos  d, 

sin  0|  = - - - '  (33) 

The  geodetic  latitude  <?  ot  the  point  is  given  by 

tan  0  =  1 .003364  tan  0i  ,  (34) 

the  ephemeris  longitude  Xf  is  given  by 

\  =  M  -  ©  (35) 

in  which  0  is  tb.e  local  hour  angle  of  the  shadow  axis,  and  the  longitude  is  given  by 
X  =  X  -  1.00  2738  isT.  (36) 

t 


*  Note  the  ioomote  lollowmg  (D2I. 


Q 


4.2  ECLIPSE  SEMI  DURATION  ATCENTER  LINE  POINT 

Calculate  the  semi-duration  of  totality  3t  the  ionospheric  central  line  point  0 
and  X  of  (34i  and  (36)  from 

s  -  Lj/n  (37) 

where 

L,  =  I23  -f0  Ian/,  (38) 

in  which 

fo  =  Pi  Ifi  cos  (d ,  -  dj )  -  y,  sin  (d,  -  d: )]  (39) 

and  where 

n  =  [(c2  -n  cos  d)J  +  (- b  )3  ]'\  (40) 

4.3  SHADOW  OUTLINE  APPROXIMATION 

The  shadow  outline,  centered  at  the  point  specified  by  (34)  and  (36),  can  be 
approximated*  by  an  ellipse  whose  semi-major  axis  is  oriented  toward  the  sun 
along  azimuth  Azq.  The  semi-minor  axis  is  given  (in  kilometers)  by 

semi-minor  axis  =  L2  (Xa6)  (41) 

where  ae  .=.  the  equatorial  radius  of  the  earth  =  6378.160  km;  the  semi-major 
axis  is  given  by 

semi-major  axis  =  L,  (3fa$Vsin  ECq  **  LJ(3fa6>)/S’  (42) 

in  which  ECq  is  the  solar  elevation.  The  approximation  is  derived  in  ( D9 1 ). 

Explicit  formulae  for  the  solar  azimuth  and  elevation  are  given  by  (D90)  and 
will  be  repeated  in  Section  1 1.3. 


*  This  approximation,  avoids  the  tetiious  point  by  point,  albeit  more  precise,  outlined  trace 
method  developed  (primarily  for  illustrative  purposes)  in  Section  C  8  ol  Appendix  C,  cf 
Comrie,  L.  J..  "Some  Computational  Problems  Arising  in  Eclipses,"  M  N.R.A.S.  87, 483  (1927) 
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5.0  CALCULATION  OF  OBSERVER  COORDINATES 


Tin'  observer  is  assumed  to  be  located  at  geodetic  latitude  0,  longitude  X 
(West  taken  .1  positive)  and  beigbt  above  sea  level  h  (in  meters). 

5.1  GEOCENTRIC  COORDINATES  OF  THE  OBSERVER 
Calculate  the  geocentric  coordinates  P  sin  0'  ^nd  p  cos  0'  from 

p  sin  0'  =  (S  +  0.15678503  x  10*6h)  sin  0 
p  cos  0'=  (C  +  0.15678503  x  10' hh)  cos  0 

where 

S  =  0.9949  7418  -  0.0016  7082  cos  2< p  +  0.0000  0210  cos  40 
C  =  1.0016  7997  -  0.0016  8208  cos  20  +  0.0000  0212  cos  40 

5.2  "FUNDAMENTAL"  RECTANGULAR  COORDINATES 

OF  THE  OBSERVER 

Calculate  the  rectangular  coordinates  of  the  observer  with  respect  to  the 


fundamental  plane  from 

£  =  P  cos  0'  sin  (h)  (47) 

r,  =  P  sin  0'  cos  d  -  P  cos  0'  sin  d  cos  (h)  (48) 

J  =  P  sin  0’  sin  d  +  P  cos  0’  cos  d  cos  (h)  (49) 

where  (y.  the  local  hour  angle  of  the  sliadow  axis,  is  given  by 

@  =  p  -  \  =  p  -  X  -  1.002738AT.  (50) 

5.3  SHADOW  RADII  AT  HEIGHT  f 


Calculate  the  pcnumbral  and  umbra!  radii  at  the  l-  igh*  s  above  the  funda¬ 
mental  plane  from 


=  -  f  tan  /, 

(51) 

-  -  f  tan  /2 

(52) 

(43) 

(44) 

(451 

(46) 


1 1 


5.4  OBSERVER  COORDINATE  VARIATIONS 
Calculate  the  hourly  variation  of  £  and  tj  from 

£  =  n  p  cos  <$>'  cos  fl 


r?  =nisind-fd;;=/ij  sin  d.  (54) 

and  calculate  the  following  combinations: 

u  =  x  -  £;  v  =  y  -  t\\  m:  -  u2  +  v2  (55) 

u  =  x  -  j;  v  =  y  -  rj;  n2  -  u2  +  v2  (56) 

D  =  uu  +  v  v.  (57) 
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6.0  CALCULATION  OF  THE  TIME  OF  MAXIMUM  ECLIPSE 
AND  OF  CONTACT  TIMES 


In  this  sect  ion,  which  deals  with  calculation  of  the  maximum  eclipse  and  con¬ 
tact  times  lor  a  specified  observation  site,  an  iterative  procedure  is  used  which  is 
fully  described  in  Sections  L).l  and  D.2  of  Appendix  D  and  then  applied  in 
Sections  I). 2,  D.3  and  D.4  of  that  appendix.  In  essence,  one  selects  an  ephemeris 
time  "close  to"  the  time  of  interest,  ai.d  then  calculates  a  formulated  correction 
to  this  selected  time.  If  necessary,  this  corrected  time  is  employed  as  the  value 
"close  to"  the  time  of  interest,  the  cycle  repeated  and  the  time  of  interest  ultimately 
established. 

6.1  TIME  OF  MAXIMUM  ECLIPSE 

Select  a  time  T0  and  -  using  elements  appropriate  to  this  time  -  calculate  the 
correction. 


t  =  -  D/n*  (58) 

Iterate  where  necessary  until  a  value  has  been  settled  upon;  the  maximum 
eclipse  time,  so  determined,  is  designated  t  . 

6.2  PENUMBRAL  (FIRST  AND  FOURTH)  CONTACTS 

Select  a  time  T0  and  -  using,  elements  appropriate  to  this  time  -  calculate 
the  correction 


in  which 


sin  tA',. 


_D 

„2 


■  i  I  cos  i p..  I 


1  U  V  -  V  u  1 

L ,  L  n  J 


(59) 


(60) 


and  where  the  negative  sign  is  selected  for  first  contact  (immersion)  and  the  posi¬ 
tive  sign  is  selected  for  fourth  contact  (emersion).  Iterate  where  necessary. 


6.3  UMBRAL  (SECOND  AND  THIRD)  CONTACTS 


Employing  elements  appropriate  to  the  time  of  maximum  eclipse  computed 
following  (58),  calculate 

j  =  _  J2.  4  Jri.  |  cos  |  (61 ) 

u  n3  '  n 

in  which 

sin  ^  =  J-  U 1  =  JL  ;  (62) 

L2  L  n  J  I_2 

note  that  Lj  |cos  ^  l/n  is  the  semi-duration  Sy  as  given  in  (D23).  The  positive  sign 
is  selected  either  for  second  contact  (immersion)  in  the  case  of  a  total  (L5  <  0) 
eclipse  or  third  contact  (emersion)  for  an  annular  (L2  >0)  eclipse;  the  negative 
sign  is  selected  either  for  third  contact  in  a  total  eclipse  or  second  contact  in  an 
annular  eclipse. 
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7.0  CALCULATION  OF  POSITION  ANGLES 


Calculate  the  position  angle  Q(  of  the  ith  contact  point  on  the  solar  limb 
measured  eastwards  from  the  north  point  (i.e from  the  hour  circle  passing  through 
the  solar  center)  by  means  of 


\  (63)* 

where  we  note  that  for  i  =  1 ,4  the  algebraic  sign  of  sin  Q.  is  that  of  Uj  whereas  for 
i  =  2,3  it  is  opposite  that  of  . 

For  some  observational  purposes,  the  position  angle  V.  measured  eastwards 
from  the  vertex  (zenith  point)  of  the  solar  limb  may  prove  more  useful.  This  is 
given  by 


V,  =  Q,  -  C,  (64) 

where  the  parallactic  angle  C(  is  given  approximately  by 

tan  q  *{,/ 1,,  (65) 

in  which  the  sign  of  sin  C;  is  that  of  (C. 


At  any  time  prior  to,  during,  or  after  the  eclipse,  expression  (63)  yields  the  position  angle  Qf 
ot  the  line  joining  the  solar  and  lunar  centers  and  thus,  the  position  angles  of  the  contact 
points  at  the  contact  times. 
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8.0  CALCULATION  OF  DIFFERENTIAL  CORRECTIONS 


8.1  PRELIMINARY  COEFFICIENTS 


Calculate 


A, 

=  —  P  COS  0'  cos  (h) 

(66) 

a2 

=  -  (SC3  +  0.1  568h  x  10'6 )  sin  0  sin  (h) 

(67) 

A  3 

=  cos  0  sin  (h) 

(68) 

B, 

=  —  £  sin  cl 

(69) 

b2 

=  (SC3  +  0.1568h  x  10~6)  (cos  0  cos  d 
+  sin  0  sin  d  cos  (7?)  ) 

(70) 

b3 

=  sin  0  cos  d  -  cos  0  sin  d  cos  (h) 

(71) 

8.2  CORRECTION  TO  CALCULATED  TIME  OF  MAXIMUM  ECLIPSE 

Apply  the  differential  correction  (in  hours)  to  the  time  of  maximum  eclipse  < 

tm  (calculated  in  Section  6.1 )  in  the  form  * 

6tm  *  Pm  5X  +  ^  ^  +  rm  6h  <7-) 

in  which  6X  and  50  are  expressed  in  minutes  of  arc,  5h  is  expressed  in  meters  and 
where* 


Pm 

=  (u  A,  +v  B,  )/n3  •  sin  1' 

(73) 

<lm 

=  (u  A 2  +v  B})/n3  -sin  T 

(74) 

1  rr 

“  (u  Aj  +vB3)/ii!  ae; 

(75) 

a®  *=  638  x  10“  meters  and  sin  1'  =  0.000  2909. 

*  p_  and  q  are  written  with  sin  V  in  the  coefficients  of  5X  and  50  in  order  to  illustrate  a  con 
venlent  wav  of  coping  with  both  the  "radian"  nature  of  pm  6X  and  q_  50  in  the  differential 
expression,  on  the  one  hand,  and  the  requirement  of  having  convenient  units  for  the  measure 
of  differential  displacements.  This  is  done  by  noting  that  since  1  radian  =  fi7°17'45''  -  3437.75 
-  206,265”,  then  V  ~  1 /3438  radian  and  t"  =  1/206,265  radian.  Further,  we  note  that  to  a 
good  approximation,  sin  1'  =  1/3438  and  sin  l”  =  1/206,265.  Hence,  if  O  '  denotes  the  num¬ 
ber  of  arc  seconds  in  8  (radians)  then  8  (radians).  =  (*'7206,265  -  8"  sin  l”  =  sin  O'  ,  and  sim¬ 
ilarly,  0  (radians)  =  673438  =  8’  sin  V  =  sin  8'. 
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8.3  CORRECTIONS  TO  CALCULATED  TIMES  OF  FIRST  AND 
FOURTH  CONTACTS 

Apply  the  differential  corrections  tin  hours)  10  the  times  of  first  and  fourth 
contacts  (calculated  in  Section  0.2)  in  the  form 

fit  =  pfiX  +  qfi0  +  rfih  (76) 

in  which  6X  and  50  are  expressed  in  minutes  of  arc,  6h  is  expressed  in  meters  and 
where 


-  (u  A|  +  vB,  )/D  •  sin  1 ' 

(77) 

=  (u  Aj  +  vBj  )/D  -sin  1' 

(78) 

=  (u  A,  +  vB3  )/D  U(p. 

(79) 

8.4  CORRECTIONS  TO  CALCULATED  TIMES  OF  SECOND  AND 
THIRD  CONTACTS 

Apply  the  differential  correction  (in  hours)  to  the  times  of  second  and  third 
contacts  (calculated  in  Section  6.3)  in  the  following  way.  Calculate 


in  which 


Next  compute 


6  K  =  ps  6X  +  q,  60  +  r^  fih 

(80) 

p5  =  (u  B,  -  v  A,)/n  -  sin  1' 

(81) 

qs  =  (u  B2  -  v  A2)/n  *  sin  1' 

(82) 

rs  =  (u  B3  -  v  A3)/naa. 

(83) 

, ,  _  (K  +  6K) 
sm0u  - 

(84) 

where,  as  in  (62),  K  and  L5  are  evaluated  at  tm  and  compute  the  corrected  semi¬ 
duration  S;  from 

STU  =  L:  Icos  4>'J! n.  (85) 

Apply  (85)  to  the  maximum  eclipse  time  3S  corrected  by  (72). 
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9.0  CALCULATION  OF  ECLIPSE  MAGNITUDE  AND  OBSCURATION 


9.1  MAGNITUDE 


The  magnitude  is  ilefined  as  the  fraction  of  the  solar  diameter  covered  by 
the  lunar  disc  at  the  time  of  greatest  phase  in  units  of  the  solar  diameter.  Thus  for 
a  partial  edif  c  the  magnitude  M,  ,*  is  given  by 


M, 


L|  -  m 
L,  +  L: 


(86) 


or,  if  Lj  is  unavailable,  by 

M,  - 

2L,  -  0.5464 


(87) 


For  a  total  eclipse  (L2  <  0)  or  ar.  annular  eciipse  (L2  >  0), 


M2 


L,  -L, 
L,  +  L2 


(88) 


9.2  OBSCURATION 

Calculate  the  fraction  S'  of  the  solar  disc  obscured  by  the  moon  from 


S'  =  (s7  A  +  B  -  s  sin  C)/ir  (89) 

in  which 

cos C  =  (Lj  +  L2  -2m3)/(L*  -  L2)  0<C«tr  (90) 

cos  B  =(L|L2  +  m7  )/m  (L|  +  L2 )  0  B  <  rr  (91) 

A  =  rr  -  (B  +  C)  (92) 

s  =  (L(  -  L2)/(L,  +L2);  (93) 


and  where  S'  =  sJ  during  the  annular  phase  and  S'  =  1  during  totality 


'  We  note  that  f»*i .  as  given  in  IS5)  and/or  (87).  ts  s  useful  expression  for  the  fraction  of  the 
solar  diameter  eclipsed  at  any  time  prior  to  second  contact  and  following  third  contact. 
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10.0  CALCULATION  OF  HELIOGRAPHIC  COORDINATES 
OF  ARBITRARY  POINTS 

10.1  HELIOGRAPHIC  COORDINATES  OF  SUB  TERRESTRIAL  POINT 

The  hcliographic  latitude  B0  and  longitude  L0  of  the  sub-terrestrial  point 


(center  of  solar  disc)  are  gotten  from 

sin  B0  =  0.12620  sin  (X^,  -  SI)  (94) 

cos  B0  cos  (Lo  -  M)=-  cos  (\q  -  L2)  (9S) 

cos  B0  sin  (Lo  —  M)  -  -  0.99200  sin  (Xq  —  £2)  (96) 

in  which 

S2  =  73°  40'  +  50"2 5 1  (97) 

M  =  292°  766  +  I4°l 843971 6  (243  0000.5  -  J.D.)  (98) 

where 


Xq  .=.  the  longitude  of  the  sun  which  isgiven  either  by  (A.1 6)  or  (A. 1 8), 

t  .=.  the  time  in  years  from  1850  A.D.  to  the  date  of  observation, 

J  .D.  .=.  the  Julian  Day  number  of  the  date  and  time  of  observation. 

10.2  POSITION  ANGLE  OF  SUN'S  AXIS  OF  ROTATION  P 

Calculate  P,  the  position  angle  (measured  eastwards  from  north  point)  oft  in- 
sun’s  rotation  axis  from 

P  =  X  +  Y  (99) 

in  which 

tan  X  =  --cosXQtane  (100) 

tan  Y  =  -  0.12722  cos  (X0  -  22)  (101) 

where 

c  .=.  the  obliquity  of  the  echptic  given  by  (A8). 
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10.3  HELIOGRAPHIC  COORDINATES  OF  ARBITRARY  POINTS 


The  lithographic  longitude  1  and  latitude  B  of  any  point  at  position  angle  0 
and  linear  radial  distance  r  on  the  solar  disc  of  lineai  radius  r^and  can  be  cal¬ 
culated  trom 


sin  B  -  sin  B0  cos  p  +  cos  B0  sin  p  cos  (P  -  0)  (102) 

cos  B  sin  (L  -  L0)  =  sin  P  sin  (P  -  0)  (103) 

where  p  can  be  calculated  from 

sin  (P  +  —  Sw)=  ~  (104) 

re  °  r© 


in  which  Sq.  the  solar  semi-diameter  (in  minutes  of  arc),  is  given  by  (A23). 

Conversely,  if  L  and  B  are  given,  0  and  r  can  be  calculated  from  the  inversion 
of  (102)  in  the  form 

cos  P  =  sin  B0  sin  B  +  cos  B0  cos  (L  -  Lo)  cos  B,  (105) 

and  from  the  application  of  (103)  and  (104). 

10.4  HELIOGRAPH  1C  COORDINATES  OF  CONTACT  POINTS 


The  ith  contact  point  at  position  angle  Q,  and  radial  distance  r  =  rQ,  such 

that  p  =  tr/2  -  Sq  and  cos  p  =  Sq  sin  1'  and  sin  p  *  1 ,  has  the  heliographic  coor¬ 
dinates  B,  and  L,  given  by  the  good  approximation 

•  tv~/p  .■  rv  t.  ^  n  \  r>  /  t  A  /  \ 

mii  voqmii  i  }  >iii  n0  t  i. r  —  v,  t  b0 


cos  B(  sin  (Lj  —  L0)  ^  sin  (P  —  Q( ) 


(107) 


*  r  and  iq  are  meant  (see  Figure  E-2)  as  linear  measures  in  some  convenient  scale  on,  say,  a 
photograph. 


11.U  TOPOCENTRIC  LUNAR  AND  SOLAR  QUANTITIES 


11.1  TOPOCENTRIC  RIGHT  ASCENSION  DECLINATION  AND 
SEMI  DIAMETER  OF  MOON 

Tin'  topocentric  right  ascension  a'( ,  declination  b'f ,  and  distance  r‘(  ol  the 
moon  are  given  by  the  following  exact  relations 

rj  cos  cos  =  rt  cos  6r  cos  qc  —  P  cos  o'  cos  r^  (10S) 

r'c  cos  sin  =  rv-  cos  6t  sin  af  -  P  cos  O'  sin  <  1 0^ ) 

r'r  sin  =  r€  sin  S<  Psin  O'  (110) 

in  which  at  and  are  the  geocentric  right  ascension  and  declination,  respectively; 
t£s  is  the  local  sidereal  iime  (L.S.T.)  at  r  hours  l-.T.  as  given  by  (84)  and  rc  is  the 
geocentric  lunar  distance  given  by 

rt  =  a^. /sin  rr  (111) 

where  the  lunar  parallax  n{  comes  from  ( A45). 

The  topoccntric  lunar  semi-diameter  Sj  is  gotten  from 

sin  Sj  =  sin  S«  (112) 

r* 

where  the  sine  of  the  geocentric  lunar  semi-diameter  is  given  by  (A48).* 

11.2  TOPOCENTRIC  RIGHT  ASCENSION  AND  DECLINATION  OF  SUN** 

The  topocentric  right  ascension  a0  and  declination  5^,  of  the  sun  are  given 
by  the  following  approximate  relations 


-  q0  -  Aa0 

(113) 

=  60  -  A60 

(114) 

*  The  difference  -  Ct  is  called  the  augmentation 

*  *  The  semi  diameter  of  the  sun  is  essentially  indifferent  to  topocentric  distance  variations. 
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in  which 


|p  cosd'  sec  6©,sin  (t£s  -  a@  )]  (115) 

lpsin0'^.s  6q  -  p  cos  <>' sin  -  <*©)’  (116) 


^a©  -  w© 

^6©  ~  ff© 
where  the  solar  parallax  itQ  is  given  by  (A22). 

11.3  SOLAR  OR  LUNAR  AZIMUTH  AND  ELEVATION; 
SOLAR  LUNAR  ANGULAR  DISTANCE 


The  azimuth  Az  (measced  east  from  north)  and  the  elevation  EC  of  either 
the  sun  or  moon  can  be  calculated  from 

cos  F.C  cos  Az  =  sin  6’  cos  <t>  -  cos  6'  sin  0  cos  (r£s  -  «')  (117) 

cos  EC  sin  Az  -  -cos  6'  sin  (r£s  -  a')  (118) 

sin  L<  =  sin  6’ sin  d+ cos  6' cos  ^  cos  (  Tgs  -  a')  (119) 

in  which  (r^  -  a')  is  recognized  as  the  hour  angle  of  the  body  in  question. 

The  solar-lunar  angular  distance  (tne  great  circle  separation  of  the  centers  of 
the  sun  and  moon)  p',  is  given  by 

cospj  =  sin  6q  sin  b'f  -!•  co  Jgcos  8^  cos  (a©  -  aj)  (120) 

or  from  (D.71)by 


P; 


2m 

L,  +  L2 


So 


(121) 
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12,0  PREDICTION  OF  SUN  SPOT  IMMERSION/EMERSION 


In  this  section,  three  alterative  methods  will  be  presented  'or  deteir.iiniig 
those  times  at  which  a  particular  point  on  the  solar  sir  face  (a  sun  spr...  f*--; 
example)  is  immersed  by  the  leading  edge  and  emerges  from  the  trailing  edge  of  the 
moon.  The  first  method  employs  the  geocentric  right  ascension  and  deeimaii.'n  of 
the  solar  location  as  coordinates  of  a  fictitious  star  and  proceeds  with  a  “standard” 
stellar  occultation  calculation.  The  s*  .ond  method  is  based  or,  the  L'.  that  at 
immersion  or  emersion  the  separation  of  the  point  on  the  s..:i  a.td  lunar  renter  is 
equal  to  the  semi-diameter  of  the  moon.  The  tliird  alternative  ih?s  a  cine- 
tnatic  simulation  composed  from  a  temporal  sequence  of  stills”  and  the  conse¬ 
quent  “visual"  observation  of  r.ot  only  the  occultation  of  selected  solar  points  but 
also  the  visual  determination  of  other  local  circumstances  as  well. 

12.1  “STANDARD"  OCCULTATION 

If  the  nehographic  coordinates  of  the  spot  1^ ,  Bs  are  given,  these  must  be 
converted  i.  j  the  geocentric  coordinates  p,  s  and  by  means  of  the  relations 

cos  ps  =  cos  B0  cos  BjCos  (Ls- L0)  +  sin  B0  sin  Bj  (122) 

sin  ps  sin  (P  -0%)  =  cos  Bs  sin  (Ls-  L0)  (123) 


and 


P..S 

sin  (ps  +  p,  ,5)  =  — - 


(124) 


The  right  ascension  as  and  declination  6,  of  the  spot  arc  calculated  according 
to 


a,  58  a®  +  Pi.s  sin  d%  sec  50 

(125) 

\  **  5©  ‘‘  Pi *  cos  O'. 

(126) 

Calculate  the  rectangular  coordinates  x,  y,  z  of  the  moon  (in  units  of  the 
equatorial  radius  of  the  earth)  with  respect  to  that  fundamental  plane,  whose  axis 
is  maintained  parallel  to  the  earth  center-sun  spot  vector,  from 

x  =  r€  (cos  b(  sin  (at  -  a$)! 

(127) 

y  =  ir  [sin  6C  cos  6s  -  cos  6t  sin  5scos  («f  -  c*s)  ] 

(128) 

z  =  rc  (sin  5f  sin  +  cos  6f  cos  5S  cos  (at  -  a5)  ] 

(129) 

in  winch 


rt  =  1  /sin  vrt  ( 1 30) 

Calculate  the  ephemeris  hour  angle  ps  of  the  axis  from 

P5  =  E.S.T.  -a  (131) 

Calculate  the  hourly  variations  x,  y,/is  and  6S  of  the  Besselian  elements  x,  y, 
ps  and  6s,  respectively. 

Calculate  the  rectangular  coordinates  of  the  observer  with  respect  to  the  fun¬ 


damental  plane  from 

£  =  p  cos  0'  sin  ©s  (132) 

r?  =  p  sin  0'  cos  -  p  cos  0'  sin  6s  cos  ©s  (133) 

f  =  p  sin  0'  sin  6S  +  p  cos  0'cos  t  cos  ©s  (134) 

in  which  ©  ..  the  local  hour  angle  of  the  axis,  is  given  by 

©r  Ps  -  -  1.002738  AT  (135) 

Calculate  the  hourly  variations  of  £  and  r?  from 

£  =  p3  p  COS0'  cos  ©s  (136) 

7  -  P,  £  sin  6?  (137) 

and  calculate  the  following  combinations 

u  =  x  -  £;  v  -  y  -tj  ;  mJ  =  u2  4  v  2  (138) 

u  =  x  -  £;  v  =  y  -  rj;  n2  =  u2  +  vJ  (1 39) 

D  -  uii  +  vv  .  (140) 


26 


In  the  spirit  of  Section  6.0  make  an  “educated  guess”  as  to  the  occultation 
time  T0  and  calculate  (138),  (139),  and  (140)  at  that  time.  Next  calculate  the 
correction  r  to  T0  from 


in  which 


0.2724880 

n 


|  cos  ^  i 


sin  \J/ 


_ I _ 

0.2724880 


(141) 


(142 


the  negative  sign  is  selected  for  immersion  and  the  positive  sign  0  selected  for 
emersion.  The  occultation  time  is  then  given  by  T0  +  t  which  can  be  used  to 
recalculate  (138),  (139),  and  (140)  in  preparation  for  a  second  “run”  through 
(141)  and  (142)  if  higher  precision  is  required. 

The  position  angles,  measured  eastward  from  the  north  point  of  the  moon 
he.,  from  the  hour  circle  passing  through  the  lunar  center),  of  immersion  and 
er.'  "ion  on  the  lunar  limb  are  given  by 

Ps  -  N  +  t 

where 

tan  N  =  u/v 

in  which  the  sign  of  sin  N  is  that  of  u. 

12.2  SUN  SPOT-LUNAR  CENTER  SEPARATION  METHOD 

Given  a  spot  having  a  position  angle  9%  and  a  linear  radial  distance  from  the 
solar  center  of  r,.,  the  spot’s  distance  from  the  center  of  the  moon  <R  is  given  at  an 
arbitrary  time  following  first  contact  and  prior  to  fourth  contact  by 

=  IP?*  -i  b3  -  2p,  s  b  cos(6s  -  Qc ) ] Vi  (145) 

p..*=  (4^)s© 
b  -  (  w)  s° 


(146) 

(147) 


(143) 

(144) 


in  which 


and  where  Qt ,  the  position  angle  of  the  lunar  center,  is  given,  as  in  (63),  by 


tan  Q4  =u /v  (148) 

The  condition  for  immersion  or  for  emersion  is  given  by  the  condition  that  the 
spot -lunar  separation  4?  is  equal  to  the  lunar  radius  rt,  or  that 


Hence,  in  order  to  determine  the  times  at  which  the  spot  is  occulted  and  then 
reappears,  ft  is  calculated  at  appropiiate  intervals  over  the  time  span  of  interest, 
followed  by  inverse  interpolation  to  these  times  when  the  condition  expressed  in 
(149)  obtains. 

Along  these  same  lines,  it  is  interesting  to  note  that  the  heliographic  longitude 
and  latitude  of  that  portion  of  the  lunar  limb  outlined  on  the  solar  disc  can  be 
readily  calculated  at  any  time  as  fellows.  For  values  of  6  such  that  p,  «£  Sq  , 
calculate 


p,  =■  b  cos  (9  Qt)  -  [rj,  -  b5  sin2  (0  -  QC)]'A‘  .  (150) 

Tne  (p|  ,6)  pairs,  so  calculated  can  then  be  employed  via  (146),  (104),  (102)  and 
(103)  to  determine  the  corresponding  (L,  B)  pairs  of  the  lunar  limb  outline. 

12.3  CINEMATIC  SIMULATION  METHOD 


In  this  method  the  eclipse  is  cinematically  simulated  by  a  temporal  sequence 
of  “stills”  which  can  be  drawn  following  calculation  of  the  solar  and  lunar,  topo- 
centric  right  ascensions  and  declinations  (or  azimuths  and  elevations)  and  of  the 
lunar  semi-diameter  (the  solar  semi-diameter  assumed  constant).  If  topocentric 
right  ascensions  and  declinations  are  employed,  the  north  point  of  the  sun  is  easily 
identified  by  means  of  the  hour  circle  “grid.”  If  azimuths  and  elevations  are 
employed,  on  the  other  hand,  the  north  point  is  fixed  by  constructing  the  parallel 
of  declination  tangent  to  the  “northern”  limb  of  the  sun  from  the  apparent  diurnal 
motion  of  the  sun  and  by  noting  that  the  north  point  on  each  “still”  is  this  tangent 
point.  Then,  from  the  appropriately  scaled  value  of  rs  (where  the  scale  is  deter¬ 
mined  by  the  selected  right  ascension/declination  or  azimuth/elevation  grid)  and 


the  value  of  8S, 
tion  visualized. 


the  spot  citii  be  drawn  on  each  “still’s”  solai  disc  and  the  occulta- 


13.0  MISCELLANEOUS  ALGORITHMS 


13.1  ANGLE  OF  SHADOW  CONE  SURFACE  AT  OBSERVATION  SITE 

The  angle  which  the  shadow  (penumbra!  or  umbral)  cone  surface  makes  at  an 
observation  site  is,  in  fact,  the  angle  which  a  generating  ray  of  the  cone  (connec¬ 
ting  the  observation  site  to  the  appropriate  contact  point  on  the  sun)  makes  at  the 
observation  site. 

Thus,  from  the  calculation  of  the  ith  contact  position  angle  Qt  in  (63)  and 
the  fact  that  Pj  (  ,  the  topocentric  angle  subtended  at  the  earth  by  the  contact  point 
and  the  solar  center,  is  the  solar  semi-diameter  S0,  the  topocentric  right  ascen¬ 
sion  and  declination  of  the  point  are  given,  following  (125)  and  ( 1 26),  by 

a'.  =  “©  +S0sinQisec6©  (151) 

6'  =  60+ S0  cos  Q,  (152) 

where  a'0  and  60  are  the  topocentric  right  ascension  and  declination  of  the  sun 
calculated  in  (113)  and  (114).  The  azimuth  and  elevation  of  the  point,  and  thus 
of  the  shadow  cone  surface,  are  then  readily  calculated  from  (1 1 7),  ( 1 1 8)  and  (1  19). 

13.2  CENTROID  OF  THE  UNOBSCURED  SOLAR  DISC 

Pointing  instructions  for  tracking  the  centroid  of  the  unobscured  solar  disc 
follow  directly  from  the  development  in  Section  D.9  of  Appendix  D.  There  (the 
solar  semi-diameter  being  employed  as  the  unit  of  distance),  (D.83)  gives  as  the 
angle  subtended  at  the  observation  site  by  the  centroid  and  the  solar  center 

-s;(A-_sin_Ac_os_A)  /_2m_\  s  (1S3) 

,  c  7T  ( 1  -  S')  \  L,  +  L2  /  ° 

where  the  symbols  are  those  of  Section  9.2.  Furthermore,  the  position  angle  0C  is 
given  by 

0c  =  Q*-ff  (154) 

where  Qc,  the  position  angle  of  the  solar-lunar  center  line,  is  given  by  (148).  Deri¬ 
vation  of  pointing  instructions  from  (153)and  (154)now  follows  the  discussion  of 
Section  13.1. 


29 


13.3  COORDINATES  OF  BORE  SIGHT  VECTCR/SPHEROID 
"LAVER"  INTERSECTION  POINT 

The  following  iterative  procedure  has  been  designed  to  calculate  the  latitude 
and  longitude  of  that  point  in  a  “layer”  at  height  h  L (in  meters)  above  the  earth 
spheroid,  which  is  bore-sighted  along  a  line  of  azimuth  A z  and  elevation  EC  from 
a  location  at  latitude  0,  longitude  X  and  height  h  (in  meters)  above  the  spheroid. 

Step  1  Calculate 

X  =  (p  cos  0')  a®  (155) 

Y  =  0  (156) 

Z  =  (psin*')a®  (157) 

from  (43)  and  (44)  and  calculate 

P  =  [X2  +  Y2  +ZMVl  (158) 

Step  II:  Calculate  the  nth  approximation  of  slant  range  (RL  from 

«J_n>  »  (159) 

where 

«<“>  =  [(p  + 0.15678503  x  10'6  hL)J  -  p 2  cos2  EC]’A 

-  p  sin  EC  ( 1 60) 

A(0)  =  0  (161) 

Step  111:  Calculate  the  nth  approximation  to  the  intersection 
point  coordinates  from 

X'Ln>=  X  +  (- cos  EC  cos  Az  sin  0  +  sin  EC  cos  0]  (162) 

Y[n)=  Y  +  fi(Ln)  [cos  EC  sin  Az]  (163) 

z[n)  =  7.  +  <R(Ln)[cos  EC  cos  Az  cos  0  +  sin  EC  sin  0)  (164) 

and  the  nth  approximation  to  its  geocentric  radius  from 
p[n)  =|(X(n))2  +  (Yln>)2  +(Z<n>)2!’'1  (165) 
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Step  IV:  Calculate  the  nth  approximation  to  the  geocentric  latitude  of 

the  intersection  point  0|<n)  from 


<n>-  1  j 


<n):;  tan' 


_ t 

?  |(X(Ln))1  +  (Y[n>)3]*  j 


zr 


(1  66) 


Step  V:  Calculate  the  nth  approximation  to  the  geodetic  latitude  of 

the  intersection  point  0{n)  from 


.  .  tan  <t> j  tn)  | 

<p\n)  =  tan  1  — - ^ - _ 

10.9933054  +  1.1  x  1 0 “ 9 hL  \ 


(167) 


Step  VI:  Calculate  the  nth  approximation  to  the  geocentric  radius  of 

the  sub-i  ltersection  point  from 

pjn)  =  0.99832707  +  0.00167644  cos  20<ln> 

-  0.00000352  cos  40<">  ( 1 68) 

and  calculate  the  nth  approximation  to  the  geocentric  latitude 
of  the  sub-intersection  point  0dn)  from 

<?p,n)  =  tan  (0.9933054  tan  0(Ln)!  (169) 

Step  VII:  Calculate  nth  approximation  to  layer  height  h(Ln)  from 

h<L")  =  ((po»(2  _(p<p">)r  sin3  (*[">  -  0p‘n>)]  w 

-  p(pn)  cos  (0(Ln)  -  0j,ln))  (170) 

Step  VIII :  Calculate  3f<n|  =  hL  -  hj"*.  If  3C(n)  >  e,  where  e  is 
preselected,  calculate 

,V<") 

A,n+"  = -  (171) 

sin  I  EC  +  (0  -  0'Ln))) 

and  returr*  to  Step  IS  for  (n  +  list  approximation;  if  Hln>  <  e 
proceed  to  Step  IX. 
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Step  IX 


Calculate  longitude  change  from  nth  approximation  as 


6Xln)  =  tan'1 


nn) 


x;n> 


(172) 


6X  >  0  is  in  an  eastward  direction.  Hence,  if  original  longitude 
is  west,  write  XL  =  X  (W)  —  6Xln) ;  if  east,  write  XL  =  X  (E)  +  6X(n] 


Step  X:  Final  results,  assuming  nth  approximation  satisfactory,  are: 


oL=*in,;\ 


XL  ±  5Xfn):  6?l  =  (R(l"> 


13.4  COORDINATES  OF  BORE-SIGHTED  OBJECT  OF  KNOWN  SLANT 
RmNGE 


The  following  iterative  procedure  has  been  designed  to  calculate  the  latitude, 
longitude  and  height  above  the  spheroid  of  an  object  with  azimuth  Az,  elevation  EC, 
and  slant  range  (St  observed  from  a  location  at  latitude  <f>,  longitude  X,  and  height  h 
(in  meters)  above  the  spheroid. 

Step  1:  Calculate 


X  =  (p  cos  </>')  ae  +  (R  { -  cos  EC  cos  Az  sin  <p 


+  sin  EC  cos  0] 

(173) 

Y  =  6(  [cos  EC  sin  Az] 

(174) 

Z  =  (p  sin  <t>‘ )  a^  +  (R  {cos  EC  cos  Az  cos  <t> 

+  sin  EC  sin  £>) 

(175) 

Calculate  the  longitude  change  6X  from 

6X  =  tan'1  ~ 

A 

(176) 

and,  as  in  Step  IX  of  preceding  section,  note  that  if  original 
longitude  is  west,  write  X  =  X(V')—  <5X;  if  east,  write  X  =  X(E)  +  6X. 

Step  III: 


Calculate  the  geocentric  latitude  of  object  <t>'0  from 


=  tan"1 


<  Z  I 
)|XJ  +  YJ  ]  *  ) 


(177) 
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Step  IV: 


Step  V: 


Step  VI: 


Step  VII: 


Step  VIII: 


Step  IX: 


and  calculate  the  geocentric  distance  of  object  r0  from 

r0  =  [X  2 -t  Y5  +  Z3  |Vi  (178) 

Calculate  the  nth  approximation  to  geocentric  latitude  of  sub- 
object  point  0'tn)  from 

<t>‘sin)  -  -  A<"  >>  (179) 

where  A(0)  =  0. 


Calculatc  the  nth  approximation  to  the  geodetic  latitude  of  the 
sub-object  point  from 


4>[n)  -  tan  1 


tan  <t>[<n)  "I 
.  0.9933054  J 


(180) 


Calculate  the  nth  approximation  to  the  geocentric  radius  of  the 
sub-object  point  from 


p<n)  =  0.99832707  +  0.001 67  644  cos  2<?<n) 

-0.00000352  cos  4*£n)  (181) 

Calculate  the  nth  approximation  to  the  height  h[n)  from 


Kn)  ~  ifo  -  (p(sn>*5  s;ji: (^n)- ^'<n>)r- 


-  p^nl  cos(^n)  t>'in)) 


Calculate 
A<n  +  ,)  =  sin’1 


j^jC 


sin  (b <n)  -  <P'.{n}) 


r o 


(182) 


(183) 


If  lA(n)  —  A(n+ 1  *  l  >  e,  return  to  Step  IV  for  (n  +  1  )st  approxi¬ 
mation;  if  I  A<nl  -  A,n  + 1 )  I  <e,  proceed  to  Step  IX. 

Final  results  assuming  nth  approximation  satisfactory  are: 

=<P[n)  :  X0  =*  ±  5X;  h0  =  h<n>. 
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13.5  SOLAR/LUNAR  VISIBILITY 


The  following  simple  criteria  enn  be  employed  to  determine  whether  the  earth 
either  partially  or  totally  blocks  the  view  of  the  sun  or  moon  from  a  high-altitude 
platform  (such  as  a  rocket  or  satellite).  We  will  assume  that  the  right  ascension  otp , 
declination  5p,  and  geocentric  radius  rp  of  the  platform  arc  known  functions  of 
time;  that  the  earth  spheroid  can  be  approximated  by  a  sphere  whose  radius  is  that 
of  the  actual  earth  spheroid  at  latitude  45°  (i.e.,  0.998331  a®);  that  the  earth’s 
atmosphere  has  the  effect  of  increasing  the  earth  radius  by  2%;  and  that  augmenta¬ 
tion  of  the  moon’s  senti-diameter  is  neglected.  Further,  let  a,  6,  S  and  k  represent 
the  right  ascension,  declination,  semi-diameter  and  parallax,  respectively,  of  either 
the  sun  or  the  moon. 

Calculate 

x  =  —  ccs  6p  sin  (ap  -a) 

(184) 

y  =  +  cos  5  sin  6p  -  sin  6  cos  6p  cos  (crp  -  a). 

(185) 

and  ,  /  0.998331  a*  \ 

sm  ^  )  , 

(186) 

and 

A,=  1 .02  (rrp  +0.998331  rr  +  S) 

(187) 

Aj=  1.02  (7rp  +0.998331  rr-S). 

(188) 

Then,  if 

x2  +  y2  >  sin2  A,  -  no  blocking  , 

(189) 

sin2  Aj  <  x2  +  y2  <  sin2  A]  -  partial  blocking , 

(190) 

x2  +  y2  <  sin2  A2  -  total  blocking. 

(191) 
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APPENDICES 


APPENDIX  A 


SOI  AR  AND  LUNAR  EPHEMERlULS 

A.O  INTRODUCTION 

This  appendix  presents  algorithms  by  which  geocentric  coordinates  of  bom 
the  sun  and  the  moon  can  he  calculated  given  the  date  (Julian  Day  number)  and 
ephenieris  tune  (1:1.)  of  interest.  Titese  algorithms,  in  the  case  of  the  sun,  are 
based  primarily  upon  the  “Tables”  of  Simon  Newcomb^7  *  and,  secondarily,  upon 
the  more  recently  published  “Tables”  of  Jean  Meeus  1 5  1 .  In  the  case  of  the  moon 
the  algorithms  are  based  primarily  upon  the  “Improved  Lunar  Ephemens”! 3  *  and 
secondarily  upon  Meeus. 

in  Newcomb,  the  celestial  longitude,  lath  ^de  and  distance  of  the  sun  is 
developed  to  high  accuracy  in  a  lengthy  series  of  i  .cular  and  periodic  terms;  in  the 
“Improved  Lunar  Lphemeris”  the  celestial  longitude,  latitude  and  parallax  of  the 
moon  is  developed  to  high  accuracy  -  following  the  Brown  lunar  theory  -  in  an 
even  more  extensive  series  of  secular  and  periodic  terms.  For  the  purposes  of  the 
current  work,  however,  where  the  accuracy  requirements  are  less  demanding, 
conceptually  modest,  but  numerically  significant  simplifications  have  been  super¬ 
imposed  on  these  developments.  These  are: 

a)  All  periodic  terms  with  coefficients  less  than  0"  1 0  have 
been  eliminated  from  both  the  solar  and  lunar  longitude 
developments; 

b)  All  periodic  terms  with  coefficients  less  than  0!'025  have 
been  eliminated  from  both  the  solar  and  lunar  latitude 
developments; 

A!!  periodic  terms  wiih  coefficients  less  than  2(J  units  in 
'he  8th  decimal  place  have  been  eliminated  from  the  solar 
distance  development; 

d)  All  periodic  terms  with  coefficients  less  than  0!'005  have 
been  eliminated  from  the  lunar  parallax  development;  and 
finally. 

c)  The  corrections  to  coefficients  of  periodic  terms  given  in 
Table  IV  of  1 3 !  have  been  eliminated 
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In  ul!  other  respects.*  the  developments  of  the  solar  and  lunar  tphemerides 
of  this  work  are  identical  to  those  given  in  these  fundamental  references  and  thus 
could  readily  yield  -  at  the  expense  of  increased  computational  effort  -  the  high 
accuracy  associated  with  their  presentation  in  The  America/!  tphemcris  and 
Xauiiedl  Almanac  by  restitution  of  the  terms  eliminated  in  a)  through  e). 

A.  1  PRELIMINARY  CALCULATIONS:  DATE  AND  TIME  Of  OBSERVATION 

Given  the  Julian  Day  number  of  the  date  of  observation  J.D  ,  calculate  the 
number  of  days  that  have  elapsed  since  January  0.5  E.T.,  1900,  the  fundamental 
epoch  with  Julian  Dry  number  241 5020.0.  This  is  denoted  by  d  and  is  given  by 

d  =  J.D.  -  2415020.0  (A.l) 

Next,  convert  the  Ephemeris  Time  of  the  observation  into  the  decimal 
fraction  of  a  day  by  means  of  the  relation 


t  - 


‘:.T. 

86,400 


(A.2) 


in  which  E.T.  the  ephemeris  time  of  the  observation,  is  expressed  in  seconds. 


Using  (A. 1)  and  (A.2).  calculate  the  fraction  of  a  Julian  century  of  36,525 
days  corresponding  to  the  interval  between  the  fundamental  epoch  and  the  date 
and  time  of  observation  by  means  of 


36,525 


t  =  d  +  t 


(A. 3) 


(A. 4) 


*  For  multi  component  terms,  the  retention  of  any  component  justifies  retention  of  the  • -afire, 

multi-component  term.  ] 

1 
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A.2  PRELIMINARY  CALCULATIONS:  FUNDAMENTAL  ARGUMENTS 

Calculate  the  fundamental  arguments  L  £2,  V,  J,  M,  Sn.  Tc,  L0>  2',  F  and 

D  of  Table  A-l ,  in  which 

L  .=.  geocentric  mean  longitude  of  the  moon 
£2  .=.  mean  longitude  of  the  moon's  ascending  node 

L'  .=.  geocentric  mean  longitude  of  the  sun 
V  .=.  heliocentric  mean  longitude  of  Venus* 

J  .=.  heliocentric  mean  longitude  of  Jupiter 

M  .=.  heliocentric  mean  longitude  of  Mars 

Sn  .=.  heliocentric  mean  longitude  of  Saturn 

T  .=.  heliocentric  mean  longitude  of  F.arth 

Lo  .=.  modified  geocentric  mean  longitude  of  tire  moon 

2  .=.  mean  anomaly  of  the  moon 

8'  =.  mean  anomaly  of  the  sun 

F  .=.  mean  distance  of  the  moon  from  the  ascending  node 

D  mean  elongation  of  the  moon  from  the  sun. 

In  Table  A-l  all  of  these  arguments  are  given  in  the  form  a  +  bt  +  ct2  +  dt3 ;  t  is 
given  for  the  date  and  time  of  observation  by  (A.4)  and  the  coefficients  a,  b,  c  and 
d  are  given  in  terms  of  revolutions'*  (denoted  by  a  superscript,  lower  case  r). 

A.T  PRELIMINARY  CALCULATIONS:  NUTATION  AND  THE  OBLIQUITY 

A.3.1  Calculation  of  Nutation  in  Longitude  and  Nutation  in  Obliquity 

Calculate  the  nutation  in  longitude  A'2  from 

=(-  1772327  -  0 " 0 !  737T)  sin  £2  + 

+  (0"20Prf  ^  0"00002  T>  sir,  2£2  t 
+  (-  1  "2709  -  C'.'OOOI  3  T)  sin(2£2  +  2F  -  2D)  + 

+  (071261  -0'.'P0031  T)  sin  8'  + 

+  (-  072037  -  0700002  T)  sin  (2£2  +  2F)  ,  (A.5) 


*  V,  3,  M,  Sn,  Te,  and  L0  refer  to  the  fixed  equinox  of  i860  In  contradistinction  to  all  Other 
arguments  which  refer  to  the  mean  equinox  of  d3te;  cf  13)  p.  288  and  |5]  p.  14  regarding 
this  point. 

”  We  note  that:  1°  =  0r0027?  77778:  V  =  0r 00004  62963,  V*  =  0r00000  07716 
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TABLE  A-1 


FUNDAMENTAL  ARGUMENTS 


a 

b 

cxIO30 

d  x  IQ3 

L 

0^75120  601080 

+  0 

r03660  11014  63356 

-  235980r 

+  1077 1 

il 

71995  354167 

- 

14  70942  28332 

+  432630 

+  1266 

\J 

.77693  521605 

+ 

273  79092  64963 

+  63044 

_ 

V 

.95019  202160 

+ 

445  03624  51095 

J 

65531  199845 

+ 

23  08089  70898 

M 

.81402  687500 

+ 

145  56470  68007 

c 

'“n 

.73852  Ml  203 

+ 

9  29437  29984 

Te 

.27499  653549 

+ 

273  78030  94025 

L 

0 

.74926  733024 

+ 

3660  09952  52418 

.82251  280093 

+ 

3629  16456  84716 

+ 1913365 

+  8203 

\ ' 

.99576  620370 

+ 

273  77785  19279 

-  31233 

-1900 

F 

.03125  246914 

+ 

3674  81956  91688 

-  663609 

-  IPO 

D 

.97427  079475 

•f 

3386  31921  S8393 

-  299023 

+  1077 
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and  calculate  the  nutation  in  obliquity  from 

Ac  =  (9"  2 1 00  +  070009 1 T)  cos  Si  + 

+  ( -CT0904  +  0!'00004  T)  cos  2S2  -t 
+  (075522  -  0!'00029  T)  cos(2S2  +  2F  -  2D)  + 

+  (0"0884  -  0" 0000 5  T)  cos  (2fi  +  2F)  (A.6) 

m  which  all  coefficients  are  given  in  arc  seconds.  It  should  be  noted  in  passing  that 
64  smaller  terms  associated  with  and  36  smaller  terms  associated  with  Ae  in 
the  complete  developments  of  these  quantities  (cf.  [2]  pp.  44  and  45)  liave  been 
eliminated. 

A.3.2  Calculation  of  Obliquity 

Calculate  the  mean  obliquity  of  the  ecliptic  eM  from 

eM  =  23°27'08"26  -  46"845T  -  070059T2  +  0"001 81T3  (A.7) 

Calculate  the  true  obliquity  of  the  ecliptic  e  from 

e  =  cK  +  Ae  (A. 8) 

A  4  SOLAR  RADIUS  VECTOR,  LONGITUDE  AND  LATITUDE 

A. 4.1  Calculation  of  Long  Period  Inequalities 

Calculate  the  long  period  inequalities  from 

5L'  =  6"40  sin  (23l°19  +  20°?T)  + 

t  f  1  .'882  -  0”016T)  sin  (57°24  +  150?27T)  + 

j  r\"  ~sr.  c  i  i°Oj  iir\o  /vt'n  , 

•  w.  -uu  Stll  (dl  .o  '  1  i  ?.  Ul  I  t 

+  0"202  sin  (31  5° 6  +  893°3T)  (A.9) 

Apply  long  period  inequalities  to  the  geocentric  mean  longitude  of  sun  L'  and 
mean  anomaly  of  sun  £'  to  yield 

L'  =  L'  +  5L'  (A. 10) 

i 

and 

K  =  r  +  6L'  (A.l  1 ) 
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A. 4. 2  Calculation  ol  the  Solar  Radius  Vector 


Calculate  the  unperturbed  solar  radius  vector  R  ftont 
R  =  1.0001  -033  0.0000  0070  T  + 

II 

+  (  0.0167  4‘)28  +  0  0000  4P9  T  +  0.0000  0126  T2 )  cos  C'c  + 

+  (  0  0001  4027  +0.0000  00  70  T)  cos  2C 

V 

+  t  0.0000  01  70  +0.0000  0001  T)  cos  3C  (A.12) 

Calculate  the  planetary  perturbations  in  the  solar  radius  vector  fiom  the 
cosine  terms  presented  in  Table  A-2.  Kach  term  is  of  the  form 

k  cos  (it.  +  j  |  Tc  -  Planet)  +  Angle) 

in  which  the  coefficient  K  is  given  in  units  of  the  8th  aecimal  place,  i  and  j  are 
tabulated  positive  or  negative  integers.  Planet  stands  for  V,  M,  J  or  Sn ,  and  Angle 
is  the  tabulated  phase  angle  in  degrees.  The  sum  of  these  planetary  terms  is 
denoted  by  ARp. 

Calculate  the  lunar  perturbations  in  the  solar  radius  vector  ARe  from 

ARC  =  3076  cos  D  +  85  cos  (D  +  C)  -  306  cos  (D  -  t)  - 

-  32  cos  (D  +  r>  +  83  cos  (D  -  S')  (A.13) 

where  the  coefficients  are  given  in  units  of  the  8th  decimal  place. 

Calculate  the  solar  radius  vector  R  from  the  expression 

R  =  Ru  (1  +  ARp+ ARt).  (A. 14) 

A.4.3  Calculation  of  the  Apparent  Solar  Longitude 

Calculate  the  equation  of  the  center  from  the  expression 

C  =(6R]0"057  -  17"240T  -  07052  T2)sin  C  + 

+  (72"33S  -  0"36l  T )  sin  2C'.  + 

+  (i:'U54  -  O'.'OOl  T)  sin  3C  (AT 5) 
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TABLE  A-2 


PLANETARY  PERTURBATIONS  IN  SOLAR  RADIUS  VECTOR 


Coefficient* 

«'c 

<T.  -  V) 

Angle 

543.2 

0 

-  1 

180° 

1575.4 

0 

-2 

0.12 

204.6 

0 

-3 

0.27 

86.6 

0 

-4 

-  0.11 

37.5 

0 

-5 

- 

21.6 

+ 1 

-  1 

175.9 

36.8 

+  1 

-2 

0.2 

200.1 

- 1 

~2 

168.5 

344.7 

- 1 

-3 

167.95 

45.1 

-1 

-4 

348.8 

44.7 

-2 

-3 

322  2 

21.6 

-2 

-4 

138  4 

32.5 

-2 

-5 

!Te  -  M) 

319.9 

34.5 

0 

+  1 

0?6 

473.6 

0 

+  2 

-  0.3 

38.7 

0 

+  3 

+  182.3 

34  8 

-1 

+  2 

+  40.9 

49.5 

-  1 

+  3 

+227.8 

110.1 

-  1 

+  4 

226.92 

24.6 

-  1 

+  5 

49.3 

24.2 

-2 

+  4 

277.6 

20.5 

-2 

+  5 

95  7 

32.0 

-2 

+  6 

(Te  -  J) 

94.9 

1627.3 

0 

+  1 

1°10 

927.0 

0 

+  1 

180.22 

64  7 

0 

+  3 

175.9 

47.9 

+  1 

+  1 

23.6 

23.7 

+  1 

+  2 

173.6 

56.2 

-  1 

+  1 

250.2 

336.0 

-  1 

+  2 

202  58 

184.9 

-  1 

+  3 

87.23 

40.1 

-2 

+  3 

103.3 

26.0 

-2 

+  4 

<Te-V 

353.9 

98.8 

0 

+  1 

0*36 

37.3 

0 

+  2 

180.1 

25  8 

-  1 

+  2 

182.6 

'In  units  of  the  8th  decimal  place. 
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Calculate  the  planetary  perturbations  in  the  solar  longitude  from  the  cosine 
terms  presented  in  Table  A- 3  in  which  the  coefficient  of  each  term  is  given  in  arc 
seconds;  the  sum  of  these  oo. sine  terms  is  denoted  by  ALp. 

Calculate  the  lunar  perturbations  in  the  solar  longitude  AL^  from 

ALV  =  67454  sin  D  +  0"177  sin  (D  +  C)  -  07424  sin  (D  -  C) 

+  O.'l  72  sin  (D  —  S'.)  (A.16) 

Calculate  the  solar  longitude  referred  to  the  mean  equinox  of  date  L'm  from 
1^,  =  L;.  +C  +  ALp  +  ALS  (A  17) 

Calculate  the  apparent  solar  longitude  (referred  to  the  true  equinox  of  date 
and  corrected  for  aberrr'.ion)  from 

X'  =  +  A  i//  -  20"496/R.  (A.  18) 

A.4.4  Calculation  of  the  Apparent  Solar  Latitude 

Calculate  the  planetary  perturbations  in  the  solar  latitude  from  the  cosine 
terms  presented  in  Table  A-4  in  which  the  coefficient  of  each  term  is  given  in  arc 
seconds;  the  sum  of  these  terms  is  denoted  by  A/Tp. 

Calculate  the  lunar  perturbations  in  the  solar  latitude  A/?;,  from 

A P's  =  07576  sin  F  -  0"047  sin  (F  -  C).  (A.  19) 

Calculate  the  apparent  solar  latitude  from  the  following  expression* 

0'=A0'p  + A0j.  (A.  20) 


*  The  expression  (A.  131  yields  the  latitude  for  the  mean  ecliptic  of  date  directly;  since  the 
aberrative  correction  is  negligible  and  the  latitude  is  unaffected  by  nutation,  the  expression's 
also  the  apparent  latitude 
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TABLE  A-3 


PLANETARY  PERTURBATIONS  IN  SOLAR  LONGITUDE 


Coefficient 

S’ 

IT,  -  V> 

Angler 

4'£58 

0 

-  1 

270° 

5.526 

0 

-2 

90.12 

0.666 

0 

-3 

90.41 

0.210 

0 

-4 

89.8 

0.116 

+  1 

-2 

90.7 

2.497 

-  1 

-2 

257.75 

1.559 

-  1 

-3 

257.96 

0.144 

-  1 

-4 

79.0 

1.024 

-2 

-3 

230.85 

0.152 

-2 

-  4 

227.4 

0.123 

-2 

-5 

49  8 

0.154 

-3 

-5 

<T#  -  M> 

214  1 

0'.‘273 

0 

+  1 

90°6 

2.043 

0 

+  2 

89  76 

0.129 

0 

+  3 

273. C 

1.770 

-  1 

+  2 

306  27 

0.425 

-  1 

+  3 

317.70 

0.500 

-  1 

+  4 

316.94 

0.585 

-2 

4  4 

185.82 

0.204 

-2 

+  5 

185  5 

0.154 

-2 

T  6 

185.0 

0.101 

-3 

4  6 

53.9 

0.106 

-3 

4  7 

IT§  -  J) 

53.3 

7M208 

0 

4  1 

gi;09 

2.731 

0 

4  2 

270  25 

0.:64 

0 

4  3 

265.2 

0.163 

4 1 

4  1 

110  2 

2.600 

- 1 

4  1 

174.77 

1.610 

- 1 

4  2 

292.60 

0556 

- 1 

4  3 

177.31 

0.210 

-2 

4  3 

<Te-V 

193.2 

0.419 

0 

4  1 

90°34 

0  198  » 

0 

+  2 

270.1 

0.320 

-  1 

4  1 

259.22 

0.112 

-  1 

4  2 

273.1 

TABLE  A-4 

PLANETARY  PERTURBATIONS  IN  SOLAR  LATITUDE 


Coefficient 

S' 

c 

<9 

1 

< 

Angli 

0:b29 

4  1 

-  1 

296° 

0.092 

-  1 

-  1 

244  6 

0.067 

-  1 

-  2 

244.8 

0.210 

-  1 

-3 

244.5 

0.031 

-  1 

-4 

65  4 

O.’l  66 

-  1 

<Te  -  J) 

4  2 

268  J6 
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A5  SOLAR  RIGHT  ASCENSION,  OECLINATION  PARALLAX 
AND  SEMI  DIAMETER 

A. 5.1  Calculation  ol  Solai  Right  Ascension  and  Declination 

Calculate  the  solar  right  ascension  o^atul  declination  60  from 
COS  6q  COS  Ci0  -  cos  X' 

cos£0sinu0  =sinX'cos<  -  14, 2d  P'  x  10'1  (A. 21) 

sin  §0  =  sin  X'  sin  r  +  44.48  (3'  x  1  O' ' . 

A. 5. 2  Calculation  of  Solar  Parallax  and  Semi-Diameter 

Calculate  the  ;olar  parallax  «q  from 

ire  =  ir0/R  =  8”794/R.  (A.22) 

Calculate  the  solar  semi-diameter  appropriate  for  eclipse  calculations*  from 

S0  =  S0/R  =  15'59"63/R  (A.23) 

A. 6  LUNAR  LONGITUDE,  LATITUDE  AND  PARALLAX 

A.6.1  Calculation  of  Additive  Terms 

Calculate  the  additive  terms**  of  Table  A-5;  these  are  designated  as: 


61. 

.=.  the  sum  of  the  8  sine  terms  in  L, 

6w 

.=.  the  sum  of  the  6  sine  terms  in  w, 

6  U 

.=.  the  sum  of  the  5  sine  terms  in  £2, 

6TC  =  6C 

.=.  the  sum  of  the  4  sine  terms  in  Te ,  0 

6  J 

.=.  the  single  sine  term  of  J, 

6Sn 

.=.  the  single  sine  term  of  , 

5tc 

.=.  the  sum  of  the  3  cosine  terms  in  y. 

*  For  other  than  eclipse  calculations,  the  adopted  value  of  the  semi  diameter  at  ur.it  distance 
So  is  16'01”18. 

**  Each  term  is  o«  the  form  K  ,'a  +  bt  +  ct1!  with  its  coefficient  K  listed  both  in  arc  seconds 
and  revolutions;  when  applied  to  the  fundamental  arguments  of  Table  A-1 ,  the  latter  unit  is 
the  more  useful. 

**  *  These  are.  in  fact,  the  long  period  inequalities  of  I  A.9). 


4b 


TABLE  A -5 


Ada' 

Coeff 

to: 

Semi  No. 

L 

1628 

-t  <5'.fe4 

l 

1629 

+  0.31 

L 

1836 

4  14.27 

L 

1638 

+  7.261 

L 

1639 

4  0.282 

L 

1646 

f  0.237 

L 

1646 

+  0.108 

L 

1618 

+  0.126 

w 

1631 

-  2.1 0 

c3 

1663 

-  0H8 

CJ 

1664 

-  2.076 

a' 

1666 

-  0.840 

CJ 

1666 

-  0.10 

eo 

1667 

—  0  533 

n 

1632 

+  0'63 

12 

1669 

+  0.17 

n 

1670 

+  95.96 

Si 

1671 

4  15.58 

Si 

1672 

+  1.86 

Te.  S' 

1633 

-  6  40 

Te.  S' 

1673 

-  0.27 

Te,  S' 

1674 

-  189 

Te.  S' 

1675 

4  0.20 

m  4  O 

7 

1676 

—  H.JiV 

7 

1677 

-  0.698 

7 

1678 

-  0.083 

J 

1634 

4  0°33 

s„ 

1635 

-  O'KZ 

ADDITIVE  TERMS 


ciont 


x  10u 

8 

64  3148' 

O'. 1422  222? 

23  9197 

.2336  3774 

1101  0802 

,5373  3431 

560  2623 

.7199  5354 

21  7592 

.4839  8132 

18  2870 

8453  6324 

8  3333 

.4035  3088 

9  7222 

,6t»b4  4$93 

167  0370' 

.1422  2222 

9  1049 

.5373  3431 

160  1851 

.7199  5354 

64  8148 

.4839  813? 

7  7160 

.5875  0000 

45  7562 

.8453  6324 

48  6111’  .1422  2222 

13  1  172  .6373  3431 


7404  3210 
1202  1605 
1435185 

7199  5354 
.4839  8132 
.5245  3688 

493  8271' 

20  8333 

145  8333 

15  4321 

.1422  2222 
,5875  0000 
.8453  6324 
.6104  3085 

333  1 790' 

53  8580 

6  4043 

.7199  5354 
4833  8132 
.5245  3688 

9  1666  6667' 

.3729  1r07 

23  0555  5556' 

.3729  1667 

b  x  10*2 

c  x  10* 

<  153  6238' 

+  123  2723 

4  191' 

-  1010  4982 

4  191 

-  1  4709  4228 

4  43 

-  1  4726  9147 

4  43 

-  1145  9387 

-  2148  8317 

-  7864  5335 

+  153  6238' 

-  1010  4982 

4  191 

-  l  4709  4228 

4  43 

-  1  4726  9147 

4  43 

+  905  01 18 

-  1145  9387 

+  153  672S' 

_  1010  4982 

4  19’ 

-  1  <709  4228 

4  43 

-  1  4726  9147 

4  43 

-  1  4716  2675 

4  43 

4  153  6238’ 

4  905  01 18 

-  .  1 45  9  387 

-  6771  8733 

_  1  4709  4228’ 

4  43 

-  1  4726  9147 

t  43 

_  1  4716  2675 

4  43 

+  292  7979' 


<  292  7979’ 


4’’ 


Apply  these  terms  tc  the  fundamental  arguments  as  follows: 


L. 

l 

=  L  +  5L 

(A. 24) 

s; 

=  C  +  6C,  Te  =1,  +  6X 

(A.25) 

j 

c 

—  J  +  5J 

(A. 26) 

S„c 

=  Sn  +  6S 

n  n 

(A.  27) 

c 

c 

-  £  i  oL- 5« 

(A. 28) 

F 

c 

=  F  +  5L  -  612 

(A. 29) 

=  D  +  6L 

(A. 30) 

A.6.2  Calculation  of  the  Apparent  Lunar  Longitude 

Calculate  the  Code  0:  Solar  Terms  in  Longitude  from  Table  A-6.  Each  term 
is  of  the  form 

Kq  sin  [i  £c  +  jCj  +  kf.  +  m  Dc  ] 

in  which  the  coefficient  lC  in  arc  seconds  is  given  in  column  1 ,  the  multiples  of  the 
fundamental  arguments  |i|  <  6;j  <  4;  |kl  <  5  and  |m!  <  8  are  given  in  columns  3 
through  6,  respectively,  and 

q  -  (1  +  2.208  x  lO'6)111  (1  -  6.832  x  10  8t)UI 

(1  +  2.708  x  10~6  +  139.978  5  7c),k|  (A.31) 

where  67  (expressed  in  revolutions)  comes  from  Section  6.i .  The  sum  of  the  117 
sine  terms  will  be  denoted  by  AL.. 

Calculate  the  Code  0:  Planetary  Terms  in  Longitude  from  the  sine  terms 
presented  in  the  latter  section  of  Table  A-6.  The  sum  of  the  26  sine  terms  will  be 
denoted  by  ALp, 

Calculate  the  lunar  longitude  referred  to  the  mean  equinox  of  date  from 

Ly  =  Ll  +  ALS+  ALp  .  (A. 32) 

Calculate  the  apparent  lur.or  longitude  (referred  to  the  true  equinox  of  date)* 

from 

Xa  =Lm  +A <1/  -  0"  1 89  sin  Cl  +0"168  sin  D.  (A.33) 

*  To  within  the  accuracy  of  this  tp’iemeris,  (A.33)  also  includes  aberration;  tor  higher  accuracy, 
specific  aberrative  corrections  listed  in  (2]  p.  109  are  required. 


TAW.E  A-A 


CODE  0  SOLAR  TERMS  IK  LONGITUDE 


S«r.  No. 

t 

c 

r 

e 

fc 

p. 

1 

+  0.127 

0 

0 

0 

+  6 

2 

+  13  902 

0 

0 

0 

+  4 

3* 

+  2369.912 

0 

c 

0 

+  2 

5 

+  1 ,9/9 

+1 

0 

0 

+  4 

6* 

+  191.953 

+1 

0 

0 

+  2 

7 

+22639.500 

+  1 

0 

0 

0 

8* 

-  4586.465 

+1 

0 

0 

-2 

9* 

-  38.428 

+1 

0 

0 

-4 

10 

0.393 

+1 

0 

0 

-u 

13 

0.289 

0 

+1 

0 

+4 

14* 

-  24.420 

0 

+1 

0 

+2 

-.5* 

-  668.146 

0 

+1 

0 

0 

16* 

-  165.145 

0 

+  1 

0 

~~2 

17 

1.877 

c 

+1 

0 

-4 

20 

+  0403 

0 

0 

0 

+3 

21* 

-  125.154 

0 

0 

0 

-1 

23 

+  0.213 

+2 

0 

0 

+4 

24 

+  14.387 

+2 

0 

0 

+  2 

25* 

+  769.016 

+2 

0 

0 

0 

26* 

-  211656 

+2 

0 

0 

-2 

27 

-  30.773 

+2 

0 

0 

-4 

28 

0.570 

+2 

0 

0 

-6 

31 

2.921 

+  1 

+1 

0 

+2 

32* 

-  103.673 

+1 

0 

0 

33* 

-  205.962 

+1 

+1 

0 

"2 

34 

4.391 

+  1 

+1 

0 

-4 

38 

+  0.283 

+  i 

-1 

0 

+4 

39 

+  14577 

+1 

-1 

0 

+2 

40* 

147.68/ 

+1 

-1 

c 

0 

6' 

•  23  475 

+1 

-1 

0 

-2 

42 

0.636 

+  1 

-1 

w 

—a 

45 

0.189 

0 

+2 

0 

+2 

46 

7.486 

0 

+2 

0 

0 

47 

8.096 

0 

+2 

0 

2 

48 

0.151 

0 

+2 

0 

-4 

52 

5.741 

0 

0 

+2 

+2 

53* 

-  411.608 

0 

0 

+2 

0 

54 

-  55.173 

0 

9 

+2 

“2 

58 

8.466 

+  1 

0 

0 

+  1 

5S 

+  18.609 

+1 

0 

0 

-1 

60 

+  3.215 

+  1 

0 

0 

--3 

63 

+  0.150 

0 

+1 

0 

+3 

84 

+  18.023 

0 

+  1 

0 

+1 

65 

+  0560 

0 

+  1 

c 

-1 
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AS  (Cont.l 


No. 

69 

70 

71 

72 

73 

76 

77 

78 

79 

83 

84 

85 

86 

89 

90 

91 

94 

95 

96 

98 

99 
102 

103 

104 

105 
108 
109* 
110 
tli 

115 

116 
118 
120 

123 

124 

125 

129 

130 

131 

134 

135 

136 

143 

144 

145 
149 


Coe«. 


*  1.060 
♦  35. 1 24 
-13.193 

-  1.187 

-  0.293 

-  0.290 

-  7.649 

-  8.627 

-  2.740 
+  1.181 
+  9.703 

-  2  4S4 
+  0  360 

-  1.167 

-  7.412 

-  Q.311 
0.757 
2,580 

+  2.533 

-  0.103 

-  0  344 

-  0992 
-45  099 

-  0.179 

-  0.301 

-  6.382 
+39.628 
+  9.366 
♦  0.202 
+  0.415 

-  2.152 

-  1.440 
+  0.384 

-  0.586 
+  1.750 
+  1.225 
+  1.267 
+  0.137 
+  0.233 

-  0.122 

-  1.089 

-  0.276 
♦  0.255 
+  0.534 
+  0  254 
+  1.938 


+  3 
+3 
>3 
+3 
+3 
+2 
+2 
+2 
+2 
+2 
♦2 
+2 
+2 
+  1 
+  1 
+  1 
+1 
+1 
+1 
0 
0 
+1 
+1 
+1 
+1 
+1 
+1 
+1 
+1 
0 
0 
0 
0 
+2 
+2 
+2 
+1 
+1 
+1 
+1 
+1 
+1 
0 
0 
0 
+4 


0 

0 

0 

0 

0 

*1 

+1 

+ 


+2 

+2 

+2 

-2 

-3 

-2 

+3 

+3 

0 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

+2 

+2 


♦  2 
0 
-2 
-4 
-6 
+2 
0 
-2 
-4 
+2 
0 
-2 
-4 
0 
-2 
-4 
+2 
0 
-2 
0 
-2 
+2 
0 


0 

+2 

-2 

0 

+2 

-4 

0 

-2 

+2 

0 

-2 

0 

0 

-2 

-2 

0 

-2 

-4 

+1 

+2 

0 

+1 

+2 

-2 

+1 

-2 

+2 

+1 

-2 

-2 

0 

0 

+1 

0 

0 

-1 

0 

0 

_o 

+1 

0 

+1 

+1 

0 

-1 

0 

+1 

-1 

0 

+1 

-1 

0 

-1 

-1 

0 

-3 

0 

+2 

+1 

0 

+2 

-1 

0 

+2 

-3 

0 

0 

0 
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TABLE  A-«  (Conti 


Sor.  No. 

Cctfl." 

V 

F 

D 

c 

c 

< 

c 

150 

-  0  952 

+  4 

0 

0 

-2 

155 

-  0.551 

»3 

-n 

0 

0 

156 

-  0.482 

•t  3 

+  i 

0 

-2 

157 

-  0  .100 

+  3 

•ti 

0 

—4 

162 

+  0.681 

+  3 

-i 

0 

0 

163 

-  0.183 

+  3 

-i 

0 

-2 

167 

-  0.297 

■*2 

+2 

0 

-2 

168 

-  0.161 

+2 

+2 

0 

—4 

172 

+  0.197 

+2 

-2 

0 

0 

173 

♦  0.254 

+2 

-2 

0 

-2 

177 

-  0.250 

+  1 

+3 

0 

-2 

186 

-  0.123 

<2 

0 

+2 

+2 

187 

-  3.99S 

+2 

0 

+2 

0 

188 

+  0.057 

+2 

0 

+2 

-2 

192 

-  0.459 

+2 

0 

-2 

+2 

193 

-  1.298 

+2 

0 

-2 

0 

194 

+  0.538 

+2 

0 

-2 

-2 

195 

+  0.173 

+2 

0 

-2 

-4 

198 

+  0.263 

n 

+1 

+2 

0 

205 

+  0.426 

+i 

+1 

-2 

-2 

209 

-  0.304 

+i 

-1 

+2 

0 

213 

-  0.372 

+i 

-1 

-2 

+2 

224 

+  0.418 

0 

0 

+4 

0 

227 

+  0.130 

+3 

0 

0 

-1 

237 

-  0.352 

+2 

-1 

0 

-1 

264 

+  0.113 

+5 

0 

0 

0 

295 

-  0.330 

+3 

0 

+2 

0 

CODE  0:  PLANETARY  TERMS  IN  LONGITUDE 

Sor.  No. 

Coeff." 

t 

Angle 

T.c 

Plenet 

< 

986 

0.822 

0 

0 

0C 

+1 

-IV 

987 

0.307 

0 

0 

179  8 

+2 

-2V 

1001 

0.348 

0 

0 

272  9 

+3 

-2V 

1002 

0.176 

0 

0 

271.7 

+4 

-3V 

1021 

0.129 

+  1 

0 

180 

-1 

+  1V 

1022 

0152 

+1 

0 

0 

+  1 

-IV 

1024 

0  127 

*1 

0 

180 

+3 

-3  V 

1061 

0.136 

0 

+2 

179.5 

+2 

-2V 

1097 

0.662 

-1 

+2 

180 

-3 

+  3V 

1099 

0.137 

-1 

+2 

0 

-2 

‘2  V 

1101 

0.133 

-1 

+2 

0 

+  1 

-IV 

11C2 

0157 

-1 

+2 

179.6 

+  2 

-2V 

5! 


TAPt.(  A -6  (Cent.) 


Sar.  No. 

Co«M ." 

S 

« 

An#* 

1172 

0643 

0 

0 

we^B 

-1 

1173 

0  187 

0 

0 

362  6 

-2 

-1 

1178 

0.166 

0 

0 

241.5 

1187 

0  144 

+1 

0 

1.0 

♦  1 

1188 

0  158 

+1 

0 

179.0 

-  i 

1189 

0  190 

+  1 

0 

180 

-2 

-1 

1208 

0167 

0 

+2 

178  5 

1225 

1.137 

-1 

♦2 

180.3 

+2 

~1 

+2 

-7 

1227 

0.211 

-1 

♦2 

178.4 

1236 

1269 

0436 

0240 

-1 

+2 

+2 

-2 

7.5 

1799 

1270 

0284 

+2 

-2 

172.5 

-2 

1279 

0.196 

0 

0 

180°2 

-2 

-1 

1283 

0.327 

0 

0 

224.4 

Wanat£ 

<1J 
+  2J 
+2J 
-1J 
+  13 
+23 

♦  1J 
-2J 
+  1  ! 
-33 
+2J 

♦  3-' 

+2M 

+2M 


*  m  Tables  A -6,  A-7  and  A  -8.  the  coefficients  of  those  terms  whose  serial  numbers 
carry  asterisks  are  subject  to  small  corrections  9™«n  in  Tab1*  ,v  of  l3)  ■  ,hK* 
corrections  have  teen  eliminated  in  this  report. 
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A, 6. 3  Calculation  of  the  Apjia'ent  Lunar  Latitude 


Calculate  the  (  <‘J,  /  .SWur  1  (  nm  in  l  oiiiiuli .  S  from  I  able  A  7  w  lu-rc  each 
tv, ui  has  the  toi  nr 

k >i  mii  lit  ■*  j  v;  +  kb  +  m  l\  ] 

aiul  ajniiti  v)  is  £i\cn  b\  (  A  311  The  sum  of  the  17b  s/nr  terms  will  he  denoted 
by  A/r 

Calculate 


S  =■  I  -r  A0s  1  \  A  4  > 

aiul  calculate  sin  S.  sin  3S  aiul  sin  5S. 

Next,  calculate 

(sin  S|  =  V {  sin  S  t  A  3M 

| sin  3S|  =  t h(  V  sin  3S  (A. 3<0 

aiul 

|  sm  5S|  =  in,  )'  sin  5S  ( A. 3  ’  I 

where 

T1,  =  <1  +  7. 708  x  10  *  -r  13‘W8  6  )  t\  3S1 

Calculate  the  Code  3  Solar  linns  in  I miitnlc.  7i C  front  Table  A-7  whcic 
each  term  Inis  the  form 

kq  cos  |  A  +  jV'  +  kl  (  +  mO(  1 . 

The  sum  of  the  hO  cosine  terms  will  be  denoted  by  7>  C. 

Next  calculate* 


A  =  1S.M‘>."700 

4  7  1  l  • 

(  A  3‘M 

11  =  33  x 

10  :A 

(  \  .401 

c  -  -t  :i«  >.  10 

J  A 

(A.4|  1 

n  =  5  3‘>d<,  x  10 

-  A. 

(  \  A  ; ) 

*  We  rvole  in  liaising  thal  A  ■),«■),  C. 

1;  1  »  1 

B  —  A.  C  —  A  a  nil  D  -  A  when*  •> , 

1 ]  1\  7  1 

■)  .  anil  1  1  ill 

the  values  ol  the  coefficients  given  in  Table  A  7  Code  6  Sofjr  Terns  m  i j(ih/dt\  Pnnotud 
Terms  having  Seiial  Numbeis  797.  798  anil  799.  lespectively 


TABLE  A-7 


COOE  1:  SOLAR  TERMS  IN  LATITUDE.  S 


Sw.  No. 

Co*ff . 

*c 

c 

397 

-  11279 

0 

0 

398 

+  2373.36 

0 

0 

399 

4  01 

0 

0 

400 

i  14. OS 

0 

0 

401 

0.13 

0 

0 

402 

+  0  60 

0 

0 

404 

+  0  25 

+  1 

0 

406 

+  6.98 

+1 

0 

407 

0.74 

+  1 

0 

408 

+  192.72 

+1 

0 

409 

-  13.51 

+1 

0 

410 

+22609.07 

+1 

0 

411 

+  3.59 

+1 

0 

412 

-  4678.13 

+1 

0 

413 

+  5.44 

+1 

0 

414 

-  38.64 

+1 

0 

415 

+  0.25 

♦1 

0 

416 

1.43 

+1 

0 

417 

0.03 

+1 

0 

418 

+  003 

+2 

0 

419 

+  1.02 

+2 

0 

420 

0.10 

+2 

0 

421 

+  14.78 

+2 

0 

422 

1.20 

+2 

0 

423 

+  767.96 

+2 

0 

424 

+  2.01 

+2 

0 

425 

-  152  53 

+2 

0 

<  091 

-  34.07 

+2 

0 

426 

427 

+7 

0 

428 

+  o.;. 

t 

429 

1.40 

+2 

0 

4?0 

0.07 

+2 

0 

431 

+  0.16 

+3 

0 

432 

+  2.96 

+3 

0 

433 

0.09 

+3 

0 

434 

+  50.64 

+3 

0 

435 

+  0.19 

+3 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


+1 

+2 

+3 

+4 

+5 

+6 

+6 

+4 

+3 

+2 

+1 

0 

-1 

-2 

-3 

-4 

-6 

-6 

—8 

+6 

+4 

+3 

+2 

+1 

0 

-1 

-2 

-3 

-4 

-5 

-6 

-e 

+4 

+2 
+  1 
0 
-1 
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TABLE  A-7  iCom.) 


5«r.  No. 

Corfi 

i 

c 

S' 

c 

Fc 

D 

c 

436 

-  16  40 

+3 

0 

0 

-2 

437 

+  005 

+  3 

0 

0 

-3 

438 

-  0  74 

■*3 

0 

0 

-4 

439 

+  003 

+  3 

0 

0 

-5 

440 

-  031 

+  3 

0 

0 

-6 

443 

*  0  30 

■♦4 

0 

0 

42 

444 

+  3  60 

+  4 

0 

0 

0 

445 

-  1.58 

44 

0 

0 

-2 

447 

-  0.03 

■•4 

0 

0 

-6 

448 

+  0.04 

+5 

0 

0 

42 

449 

+  0.28 

+5 

0 

0 

0 

450 

-  0.14 

+5 

0 

0 

-2 

452 

-  0.06 

0 

41 

0 

46 

454 

-  1.59 

0 

41 

0 

44 

465 

+  0.53 

0 

4  1 

0 

43 

456 

-  25.10 

0 

4  1 

0 

42 

457 

+  17.93 

0 

+1 

0 

4  1 

458 

-12698 

0 

4  1 

0 

0 

459 

4  0  32 

0 

4  1 

0 

-1 

460 

-lto.06 

0 

41 

0 

-2 

46t 

4  0.29 

0 

*1 

0 

-3 

462 

-  6.46 

0 

4 1 

0 

-4 

463 

-  0.22 

0 

41 

0 

-6 

464 

-  0.04 

0 

42 

0 

+4 

465 

-  1.68 

0 

42 

0 

42 

466 

-  0.04 

0 

42 

0 

41 

467 

-  066 

0 

+2 

c 

0 

468 

-  004 

0 

42 

0 

-1 

469 

-  16.35 

0 

42 

0 

-2 

471 

-  065 

0 

42 

0 

-4 

472 

-  0.57 

0 

' " 

0 

-2 

475 

-  0.50 

41 

4) 

0 

44 

476 

+  0.08 

4  1 

4  1 

0 

43 

477 

-  11.75 

+1 

4  1 

0 

+2 

478 

+  1.52 

+1 

41 

0 

+  1 

479 

-115  18 

4! 

4  1 

0 

0 

480 

-  0.12 

+  1 

4  1 

0 

-1 

481 

-182  36 

4  1 

4  1 

0 

-2 

482 

+  0  36 

+1 

4 1 

•  « 

0 

r. 

-3 

—a 

483 

-  8.66 

4l 

T  1 

V 

485 

-  037 

1 1 

4  1 

0 

-6 

486 

-  0.09 

+2 

42 

0 

0 

487 

-  0.27 

+2 

42 

0 

-2 

488 

-  0.16 

42 

42 

0 

-4 

490 

-  0.09 

-1 

4  1 

0 

46 

492 

-  2.27 

-1 

4  1 

0 

44 
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TAbwE  A  7  (Com.) 


S*r.  No. 

Cc«U 

493 

+  0  38 

494 

-  23.59 

495 

-  0  55 

49G 

-138.76 

497 

>  0.33 

498 

-  31.70 

499 

+  004 

600 

-  1.53 

601 

-  0.06 

602 

-  0.04 

603 

-  0.21 

604 

-  022 

606 

-  0  2. 

506 

-  007 

607 

-  1  45 

608 

4  0.14 

509 

-  1056 

611 

-  7.59 

612 

+  007 

513 

-  2.54 

614 

-  025 

615 

+  022 

516 

+  3  32 

517 

-  004 

618 

4  11.67 

619 

-  037 

620 

-  1.17 

621 

+  0.04 

622 

+  0.20 

623 

+  0.06 

-  0.17 

526 

-  094 

527 

-  0.67 

528 

-  0.08 

529 

-  0.06 

531 

+  036 

532 

+  096 

633 

-  0.23 

542 

-  0.13 

543 

-  1.25 

544 

-  6.12 

545 

-  0  65 

546 

-  003 

547 

-  007 

548 

-  2  40 

649 

-  2  32 

s 

S’ 

F 

D, 

c 

c 

c 

€ 

♦  1 

0 

+3 

-1 

<1 

0 

42 

-1 

+  1 

0 

+1 

-1 

♦  1 

0 

0 

-1 

+  1 

0 

-1 

-1 

+  1 

0 

-2 

-1 

4  1 

0 

-3 

-1 

+  1 

0 

~4 

-1 

+  1 

0 

-6 

-2 

+2 

0 

+4 

-2 

+2 

0 

+2 

-2 

♦  2 

0 

0 

-2 

+2 

0 

-2 

+2 

+  1 

0 

+4 

+2 

+  1 

0 

+2 

♦2 

♦  1 

0 

+  1 

♦2 

+  1 

0 

0 

♦2 

+  1 

0 

-2 

♦2 

4-1 

0 

-3 

+2 

+  1 

0 

-4 

♦2 

+  1 

0 

—6 

+2 

-1 

0 

+4 

+2 

-1 

0 

♦2 

+2 

-1 

0 

+  1 

+2 

-1 

0 

0 

♦2 

-1 

0 

-1 

+2 

-1 

0 

> 

A 

+2 

-1 

0 

-3 

+  2 

-1 

0 

-4 

+2 

-1 

0 

—6 

<3 

4’ 

0 

*2 

+3 

+  1 

0 

0 

+3 

+  1 

0 

-2 

+3 

+  1 

0 

—4 

+3 

+  1 

0 

-6 

+3 

-1 

0 

+2 

+  3 

-1 

0 

0 

+3 

-1 

0 

-2 

+  1 

+2 

0 

+2 

+  1 

+2 

0 

0 

rl 

+2 

0 

-2 

+  1 

+2 

0 

-4 

•o 

+2 

0 

-6 

-1 

+2 

0 

+4 

-1 

+2 

0 

+2 

-1 

+2 

0 

0 
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TABLE  A-7  (Coot.) 


S»f  No. 

C 

c 

f 

< 

F 

c 

Dc 

550 

-  1  82 

-1 

4 2 

0 

-2 

551 

-  0.12 

-1 

+  2 

0 

653 

-  022 

+  3 

0 

-2 

556 

-  004 

-1 

■*3 

0 

0 

638 

-  006 

-1 

■*3 

0 

-2 

660 

-  004 

0 

0 

•*2 

42 

561 

-  020 

0 

0 

<2 

0 

662 

4  084 

0 

0 

42 

-1 

663 

-  52.14 

0 

c 

■*2 

-2 

664 

•*  0.26 

0 

0 

■42 

-3 

665 

-  1  67 

0 

0 

42 

666 

-  o.o: 

0 

0 

42 

-6 

669 

+  0.07 

+  1 

0 

42 

-1 

670 

-  9.52 

+  1 

0 

42 

-2 

671 

1  0.04 

+  1 

0 

+2 

-3 

672 

-  0.33 

+1 

0 

+  2 

4 

674 

-  0.04 

+  1 

0 

42 

-6 

577 

-  c.-n 

-1 

0 

42 

42 

678 

♦  0.06 

-1 

0 

+2 

41 

579 

-  85.13 

-1 

0 

42 

0 

680 

+  004 

-1 

0 

42 

-1 

681 

+  3.37 

-1 

0 

42 

--2 

683 

+  0.04 

-1 

0 

42 

-4 

585 

-  0.75 

+2 

0 

42 

-2 

688 

-  1.14 

-2 

0 

42 

+2 

689 

-  0.74 

-2 

0 

+2 

0 

690 

+  0  38 

-2 

0 

42 

-2 

693 

-  0.04 

+3 

0 

42 

0 

594 

-  007 

+3 

0 

42 

-2 

595 

-  008 

-3 

0 

42 

44 

597 

-  Oil 

-3 

0 

42 

0 

COO 

+  0.04 

-3 

’2 

-  4. 

603 

+  0.10 

0 

+1 

42 

0 

604 

-  2.26 

0 

+1 

42 

~2 

606 

-  0.1/ 

0 

+1 

42 

~A 

607 

+  0.04 

0 

-1 

42 

*2 

608 

+  0.16 

0 

-1 

+2 

0 

609 

-  0.06 

0 

-1 

42 

-1 

610 

+  1.30 

0 

-1 

42 

-2 

611 

+  0.08 

0 

-1 

42 

-4 

612 

-  009 

0 

+2 

4? 

-2 

617 

-  035 

+  1 

♦1 

42 

-2 

618 

-  0.03 

41 

+  1 

42 

-4 

619 

-  007 

-1 

-1 

42 

42 

620 

+  0.31 

-1 

-1 

42 

0 

624 

+  003 

+  1 

-1 

42 

0 
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TABLE  A  7  (Cont.) 


Ser.  No. 

CoeH." 

*« 

S' 

c 

?c 

Dc 

625 

♦  007 

+  1 

-1 

+2 

4. 

627 

-  0  33 

-  \ 

+  1 

4  2 

0 

C29 

+  0.19  -1  •<' 

COPE  2:  SOLAR  TERMS  IN  LATITUDE. 

+2 

7,C 

-2 

641 

-  0.725 

0 

0 

0 

41 

642 

+  0  601 

0 

0 

0 

+2 

643 

+  0.394 

0 

0 

0 

+3 

646 

-  0.042 

c 

0 

0 

+6 

650 

-  0445 

-1 

0 

0 

44 

651 

4  0.068 

+1 

0 

0 

i  3 

652 

+  0.029 

+1 

0 

0 

42 

653 

+  0.455 

+1 

0 

0 

41 

654 

0.079 

+1 

0 

0 

0 

655 

-  0.P34 

+1 

0 

0 

-1 

656 

-  0.077 

+1 

0 

0 

-2 

657 

0.192 

+1 

0 

0 

-3 

660 

-  0.092 

+  1 

0 

0 

-6 

663 

-  0.074 

+2 

0 

0 

44 

666 

+  0.054 

+2 

0 

0 

4-1 

667 

+  0.107 

+2 

0 

0 

669 

+  5.679 

+2 

0 

0 

-2 

670 

-  0.030 

+2 

0 

0 

-3 

671 

-  0.308 

+2 

0 

0 

—t 

673 

-  0.074 

+2 

0 

0 

-€ 

676 

-  0.166 

+3 

0 

0 

42 

678 

-  1.300 

+3 

0 

0 

C 

680 

♦  O.XFo 

+3 

0 

0 

-2 

682 

4  0.042 

+3 

0 

0 

-i 

687 

-  0.145 

+4 

0 

0 

C 

688 

+  0.062 

+4 

0 

0 

~4 

696 

+  0.123 

0 

+1 

0 

u 

697 

-  0.032 

0 

+1 

0 

698 

+  0.040 

0 

+1 

0 

700 

-  1.302 

0 

+1 

0 

702 

+  0.054 

0 

+1 

0 

—4 

703 

+  0.031 

0 

+1 

0 

n 

J-1 

o 

/U4 

—  U.HID 

w 

707 

+  0.131 

0 

+2 

0 

+ 

708 

-  0.037 

0 

t2 

0 

709 

-  0.740 

0 

+2 

0 

7 1 1 

-  0.044 

0 

+2 

0 

712 

-  0.C25 

0 

43 

0 

“ 

715 

+  0.041 

+1 

+ 1 

0 

+ 

i 
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TABLE  A-7  iCont.) 


Ser.  No. 

CoeH." 

£ 

c 

c 

F 

c 

D 

c 

7  ’7 

+ 

0.787 

4  1 

+1 

0 

+2 

719 

+ 

0.461 

4-  1 

4  1 

0 

0 

721 

-4 

2.056 

+  1 

+  1 

0 

-2 

723 

- 

Ox-'! 

+  1 

+  1 

0 

— 

725 

- 

0.027 

+  1 

4-1 

0 

-6 

731 

+ 

0.146 

-1 

+  1 

0 

4-4 

733 

- 

0.443 

-1 

4-1 

0 

42 

735 

+ 

0.670 

-1 

+  1 

0 

0 

737 

- 

1.540 

-1 

4-1 

0 

-2 

739 

0.111 

-1 

+  1 

0 

-4 

744 

+ 

0.116 

+2 

+  1 

0 

42 

746 

+ 

0.259 

+2 

+  1 

0 

0 

747 

+ 

0.078 

+2 

4-1 

0 

-2 

752 

- 

0212 

+2 

-1 

0 

42 

753 

- 

0.151 

+2 

-1 

n 

0 

760 

+ 

0.032 

+3 

4-1 

0 

0 

766 

- 

0.026 

+3 

-1 

0 

0 

777 

+ 

0.117 

+  1 

+2 

0 

-2 

778 

- 

0.032 

+1 

4-2 

0 

-4 

782 

+ 

Q.027 

-1 

4-2 

0 

0 

783 

- 

0  105 

-1 

♦2 

0 

-2 

CODE  3:  SOLAR  TERMS  IN  LATITUDE,  N 

787* 

-526.G69 

0 

0 

+1 

-2 

788 

-  3.352 

0 

0 

+1 

-4 

789* 

4  44.297 

4 1 

0 

41 

-2 

790 

-  6.000 

4 1 

0 

41 

-4 

791* 

4  20.599 

-1 

0 

41 

0 

792 

-  30.598 

-1 

0 

4  1 

._o 

793 

-  24  649 

-2 

0 

+ 1 

0 

794 

-  2.000 

-2 

0 

41 

-2 

795 

-  22.571 

0 

4 1 

41 

-2 

79C 

4  10.985 

0 

—  1 

4 1 

-2 

CODE  6:  SOLAR  TERMS  IN  LATITUDE,  PRINCIPAL  TERMS 

7S7  +  18518'.^11  sinS 

798  +  1.189  sin  S 

799  6.241  sin  3S 
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TABLE  A-7  (Coni.) 


CODE  4;  PLANETARY  TERMS  IN  LATITUDE 


S«r.  No. 

Cocff." 

«c 

F 

4. 

Cc 

1428 

+0.045 

0 

-1 

+2 

1437 

+0.068 

0 

-1 

+2 

1442 

+0.029 

-1 

+1 

+2 

1450 

+0.031 

-1 

-1 

+2 

1468 

+0.027 

0 

0 

0 

1472 

+0.077 

0 

0 

0 

1476 

+0.025 

0 

0 

0 

1477 

+0.074 

0 

0 

0 

1483 

+0.030 

0 

0 

0 

1534 

+0.051 

+1 

+1 

-2 

1535 

+0.051 

+1 

-1 

-2 

1546 

+0.035 

0 

0 

0 

1562 

+0083 

0 

+1 

0 

Lo 

ft 

Angi« 

T.c 

Planete 

0 

0 

0° 

-3 

+3V 

0 

0 

270 

-6 

+5V 

0 

0 

180 

-3 

+3V 

0 

0 

180 

-3 

+3V 

+1 

0 

285 

-I 

+  1V 

+1 

0 

215.6 

+5 

-3V 

+1 

0 

255 

-6 

+4V 

+1 

0 

51,6 

-6 

+3V 

+1 

0 

125 

+8 

-0V 

0 

0 

0° 

-2 

+2J 

0 

0 

0 

-2 

+2J 

+1 

0 

168 

0 

+2J 

0 

+2 

0 

0 

0 
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Calculate  the  Code  .?  Solar  Terms  in  Latitude,  JV  from  Table  A-7  where  each 
term  has  the  form 

Kq  sin  li¥c  +  jfci'  +  k  f;c  +  m  Dc  ] . 

Tlie  sum  of  the  lOsrm  terms  will  be  denoted  by  N. 

Calculate  the  Code  4:  Planetary  Terms  in  Latitude  from  Table  A-7;  the  sum 
of  the  13  sine  terms  (each  similar  in  form  to  the  planetary  perturbations  of 
Section  A.4.2)  will  be  denoted  by  Af}p. 

Calculate  the  apparent  lunar  latitude  (referred  to  the  true  ecliptic  of  date)  of 
the  center  of  mass  of  the  moon  from 

(3  =  A  (sin  S]  +  B  (sin  3S)  +C  (sin  5s)  +  DN  +  A0P+  (A.43) 

+  0'::i5  sin  L 

Tor  use  in  eclipse  calculations  only,  calculate  the  apparent  latitude  of  the 
center  of  figure  of  the  moon  (3  (F)  from  the  expression 

0  (F)  =  p  -  0!’ 6  (A. 44) 

A.6.4  CALCULATION  OF  THE  LUNAR  PARALLAX 

Calculate  the  Code  .V  Solar  Terms  in  Parallax  and  the  Code  5.  Planetary 
Terms  in  Parallax  from  Table  A-8.  The  sum  of  the  70  cosine  solar  terms  and  the 
2  cosine  planetary  terms  will  be  denoted  by  sin  rr  (in  arc  seconds). 

Calculate  the  horizontal  parallax  of  the  moon  rrt  from 

itt  =  sin  rr  (0.9999  53253  +  (3.9168  x  10  (si”*)2 )  - 

-  0"0890  -  0:'0049  cos  C.  (A. 45) 

where  »rt  is  in  arc  seconds. 

A. 7  LUNAR  RIGHT  ASCENSION,  DECLINATION  AND  SEMI-DIAMETER 
A.7.1  Calculation  Lunar  of  Right  Ascension  and  Declination 

Calculate  the  lunar  right  ascension  at  and  declination  5t  from 

cos  6t  cos  oct  =  cos  P  cos  X 

cos  6,  sin  Oy  =  cos  p  sin  X  cos  e  -  sin  P  sin  c 

sin  5t  =  cos  j3  sin  X  sin  e  +  sin  p  cos  e  (A. 46) 
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TABLE  A-8 


Sv  No 

CODE  5:  SOLAR  TERMS  IN  PARALLAX 

coetf."  s; 

F 

c 

°c 

802 

+  0.7007 

0 

0 

0 

+  4 

803' 

*  28.2333 

0 

0 

0 

+  2 

804 

+  3422.7030 

0 

0 

0 

0 

808 

+  0.0433 

+  1 

0 

0 

+  4 

807 

+  3.036' 

+  1 

0 

0 

+  2 

808 

+  186.5398 

+  1 

0 

0 

0 

809" 

+  34.3117 

+  1 

0 

0 

-2 

810 

+  0.6LK>3 

-r  "l 

0 

0 

-4 

811 

+  0.0086 

+  1 

0 

0 

-6 

81 0 

-  0.0063 

n 

+ 1 

0 

+  4 

814 

-  0.3000 

0 

+ 1 

0 

+  2 

816 

-  C.3997 

0 

+ 1 

0 

0 

816 

+  1.S178 

0 

+ 1 

0 

-2 

817 

+  0.0339 

0 

+ 1 

0 

-4 

820 

--  0.9781 

0 

0 

0 

+  1 

821 

+  0.0054 

+  2 

0 

0 

+  4 

822 

+  0.2833 

+  2 

0 

0 

+  2 

823 

+  10.1657 

+  2 

0 

0 

0 

824 

-  0.3039 

+  2 

0 

0 

-2 

825 

+  0.3722 

+  2 

0 

0 

-4 

826 

+  3.0109 

+  2 

0 

0 

-6 

829 

-  0.0484 

+  1 

+  1 

0 

+  2 

830 

-  0.9490 

+  1 

+  1 

0 

0 

831 

♦  1.4437 

+  1 

+  1 

0 

-2 

832 

♦  0.0673 

+  1 

+ 1 

0 

-4 

834 

+  0.0060 

+  1 

-1 

0 

+  4 

835 

♦  0  2302 

+  1 

-1 

0 

+  2 

836 

+  1.1523 

+  1 

-1 

0 

0 

837 

-  0.2257 

+  1 

- 1 

0 

-2 

838 

-  0.0102 

+  1 

-1 

0 

-4 

841 

-  0.0085 

0 

+  2 

0 

0 

842 

+  0.0916 

0 

+  2 

0 

—  2 

845 

-  0.0124 

0 

G 

+  2 

C 

846 

-  0.1052 

0 

0 

+  2 

-2 

849 

-  0.1093 

+ 1 

0 

0 

+  1 

850 

+  0.0118 

+■ 1 

0 

0 

- 1 

851 

-  0.0386 

♦  1 

0 

0 

-3 

854 

+  0.1494 

0 

+ 1 

0 

+  1 

858 

+  0.0243 

+  3 

0 

0 

+  2 

859 

+  0.6215 

+  3 

0 

0 

0 

860 

-  0.1187 

+  3 

0 

0 

-2 

861 

+  0.0074 

+  3 

0 

0 

-4 

854 

-  0.0051 

+  2 

+ 1 

0 

+  2 

86S 

-  0.1038 

+  2 

+ 1 

0 

0 

866 

-  0.0192 

+  2 

+ 1 

0 

“2 

867 

+  0.0324 

+  2 

+ 1 

0 

-4 

870 

+  0,0213 

+  2 

-1 

0 

+  2 

871 

+  0.1268 

+  2 

-1 

0 

0 

875 

-  0.0106 

+ 1 

+  2 

0 

0 

876 

+  0.0484 

+  1 

+  2 

0 

-2 
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TABLE  A -8  (Continued) 


Ser.  No. 

Coetf.” 

f 

c 

F 

c 

D 

c 

880 

+  0.0112 

+  1 

-2 

0 

+  2 

881 

+  0.0196 

+  1 

-  2 

0 

0 

882 

-0.0212 

+  1 

-2 

0 

-2 

888 

-0.0833 

+  1 

0 

+  2 

-2 

892 

-0.0481 

+  1 

0 

-2 

+  2 

893 

-0  7136 

+  1 

0 

-2 

0 

894 

-0.0112 

+  1 

0 

-2 

-2 

896 

-0.0066 

0 

+  1 

+  2 

-2 

900 

-0.0100 

+  2 

0 

0 

+  1 

901 

+  0.0155 

+  2 

0 

0 

-  1 

902 

-0.0088 

+  2 

0 

0 

-  3 

905 

(•0.0164 

4  1 

f  1 

0 

+  1 

911 

+  0.0071 

0 

0 

+  2 

-  1 

914 

+  0.0401 

+  4 

0 

0 

0 

915 

-0.0130 

+  4 

0 

0 

_  i 

* 

918 

-0.0097 

+  3 

+ 1 

0 

0 

923 

+  0.0115 

+  3 

- 1 

0 

0 

939 

-0.0090 

+  2 

0 

+  2 

-2 

941 

-0.0053 

+  2 

0 

-2 

+  2 

943 

-0.0141  +2  0  -2 

CODE  5:  PLANETARY  TERMS  IN  PARALLAX 

-2 

Ser.  No. 

Coeff 

«.  °c 

Angle 

T.c 

Pl«tet( 

1580 

+  0.0055 

-1  +2 

180° 

-3 

+  3V 

1610 

+  0.0095 

- 1  +2 

180°. 3 

+  2 

-  2J 

63 


in  which  ('  is  giu'ii  either  by  (A. 4?)  or  (A. 44)  depending  upon  whether  ■‘usual" 
or  eclipse  calculations  are  contemplated. 

A. 7. 2  Calculation  of  Lunar  Semi  Diameter 

Calculate  the  lunar  semi-diameter  (in  arc  seconds)  from  the  expression 

St  =  0.07W  +  0. 272453  (A  47) 

or  from 

sin  Sj  =  0.272488  sin  »rt  (A.48) 


M 


APPENDIX  B 


EPHEMERIS  SIDEREAL  TIME 

11k-  I  phemeris  Sulereal  lime  (i.e.,the  right  ascension  of  llu'  I  phemeris  Min 
ili;m)  al  (l1'  I  .1  .  on  tire  il.ite  of  I  lie  observation  is  calculated  from  llu*  expression 

I  .S  r.(0h  >  -  (V.:7(>')|  ').W7<>5  ■*  O’.OOJt;  7»>1  s.l  s>7»V4  7il  r 

+  8U5.M-'  x  10  0 ■'  (III) 

whcie  il  is  given  by  (A. I).  Tlu1  increment  ol  sidereal  time  associated  with  the 
interval  extending  Irom  01’ l.T.  to  the  lime  of 'observation  is  given  b\ 

a  i  s. i .  =  ir.w:7;>  7')iv>:  o5  +  ti-i.y  x  io  "  di  r  ai.:) 
where  r  is  defined  by  (A. 2). 

The  "equation  ol  the  equinoxes"  (A^  cxn  i)  is  added  next  to  I  he  sum  of 
(B.U  and(IL2)  to  yield  Iheapparcnt  ephemeris  sidereal  time  on  the  ilale  a  nil  at  the 
(ephemeris)  time  of  observation 

I  S  1 .  =  li.ST.  (0h)  +  A  l-'.S.T.  +  0  L>  1  74  Av  t»..i> 

where  Ai \>  is  given  In  (A. 51. 

I  inallx ,  for  an  observer  at  longitiule  X  (West  ivsitive;  Iasi  negative  I  the  loeal 
siilereal  time  (I..S  T.)  at  r  hours  l  .T.  is  given  by 

LSI  =  1ST.  X  (11.4)* 


*  We  note  that  the  Ephemeris  Sulereal  Time  ol  (B.3I  at  r  horns  tphemei  is  T  mie  is  nmnet  ic.illy 
•quel  to  Greenwich  Sulereal  Time  (i  e  ,  the  right  ascension  ol  the  Greenwich  Memlian)  at  r 
hours  Universal  Time,  the  same  comment  obtains  lor  lB  41  also 
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APPENDIX  C 


PREDICTION  OF  SHADOW  OUTLINE  AND 
MOTION  ON  THE  SPHEROID 


C.O  INTRODUCTION 

This  appendix  develops  the  algorithms  necessary  to  predict  the  outline  and 
motion  of  the  moon's  shadow  on  an  earth  spheroid  of  aibitrary  radius  and  flatten¬ 
ing.  The  development  is  based  upon  t'hauvenet's  expositon  of  Bessel's  original 
treatment  ( 1  ]  and  later  summarized  authoritatively  in  “The  Explanatory 
Supplement"  1 2 1 

The  development  is  being  repeated  here  in  detail  because  the  summary  of 
requisite  formula  in  [21  with  an  ordering  quite  different  from  that  of  [11  lacks 
both  the  justification  and  seeming  inevitability  with  which  elements  related  to  the 
flattened  spheroid  and  the  sluidow  motion  thereon  were  originally  introduced. 
Furthermore,  it  seems  appropriate  to  show  (in  view  of  today’s  interest  in  upper 
atmospheric  and  ionospheric  research  and  the  employment  of  high  altitude  and 
orbiting  instrument  platforms)  that  the  relationships  pertaining  to  center  line  posi¬ 
tion  and  duration  in  the  ionosphere  -  originally  cited  by  Lewis*  in  1940  -  follow 
naturally  from  a  development  in  which  scale  factors  arc  explicitly  employed. 

C.1  SHADOW  AXIS  COORDINATES  AND  THE  FUNDAMENTAL  PLANE 

Given  the  right  ascension  o0.  declination  60  and  parallax  n0  of  the  sun  and 
the  right  ascension  at ,  declination  6t  (calculated  from  the  lunar  latitude  of  the 
center  of  figure)  and  parallax  of  the  moon,  we  will  calculate  the  right  ascension 
a  and  declination  d  of  the  solar -lunar  shadow  axis  from  the  equations. 

G  cos  d  oos  a  =  r0  cos  60  cos  a0  -  rt  cos  6t  cos  a4 

G  cos  d  sin  a  =  r0  cos  60  .dn  o;0  -  rt  cos  o4  sin 

G  sin  d  =  r0sin  6 q  —  r4  sin  (CM) 

where  G  is  the  solar-lunar  separation.  Following  division  by  r0.  (CM)  can  be 
written  in  the  form, 


g  cos  d  cos  a  =  cos  6 0  cos  a0  -  b  cos  6C  cos  at 

g  cos  d  cos  a  =  cos  60  sin  a0  -  b  cos  6r  sin  af  (CZ) 

g  sin  d  =  sin  50-  b  sin  6f 

*  Lewis,  I.  M.  "Formulas  tor  the  Ionospheric  Track  in  Eclipses,"  Astron  J  40.  4,  (1940) 


(>7 


(C.3) 


in  w  iucli  g  =  C/rw  •  is  given  by 

r  a0,/sin  rr4  sin  »Tq 

h  =  -k~  * - =  - 

•\h  a^./sin  ttq  sin  rrt 

ami  a,it  is  t tic  equatorial  radius  of  t tie  earth. 

C.2  COORDINATES  OF  THE  MOON  RELATIVE  TO  THE 
FUNDAMENTAL  PLANE 

Next,  following  Bessel,  we  introduce  a  right  handed  cartesian  system  with  its 
original  theearth’soenterand  its  7-axis  maintained  parallel  to  the  (moving)  shadow- 
axis  of  right  ascension  a  and  declination  d  as  determined  by  (C.2)  and  (C.3).  The 
xy -plane,  so  determined,  is  the  fundamental  plane;  the  x-axis,  which  is  positive  to 
the  cast,  is  traced  out  by  the  intersection  of  this  fundamental  plane  and  the  equa¬ 
torial  plane;  the  y-axis,  which  is  positive  tr  ward  the  north,  completes  the  right 
handed  triad  of  this  system  portrayed  in  Figure  C-l. 

The  cartesian  coordinates  of  the  moon  in  this  system  follow  immediately 
from  application  01 an  Euler  Angle  transformation*  to  the  lunar  coordinates  in  the 
geocentric  equatorial  system  After  some  trigonometric  reduction,  these  coordin¬ 
ates  become 

x  =  *a0 


y  =  v  a* 

z  =  ca<f, 

in  which 


(C.4) 


*  =  (cos  sin  (a*  -  a)  j  /sin  »rc 

>'  ~  (sin  6C  cos  d  -  cos  6C  sin  d  cos  (af  -  a))  /sin  »rt 

z  =  (sin  6t  sin  d  +  cos  6t  cos  d  cos  (ctt  -  a))  /sin  (C.5) 


*  In  this  transformation  we  successively  rotate  positively  (eounierckickwise)  about  the  original 
zaxis  by  <P,  positively  about  the  new  x  axis  by  0,  and  positively  about  the  new  z-axis  by  ^ ; 
the  rotation  matrix  RU',0.0)  so  obtained  relates  the  old  (x.y.z)  system  to  the  new  (x  .y'.z  ) 
system  via  the  matrix  equation 


x’ 

Y' 

z' 


R(  ip  At) 


for  the  case  at  hand  tv  *  0, 


*  -  -  d.  and  <P  =  a  +  — 
2  2 


68 


■iiid  wlr  it  should  be  noted  that  l  he  \  amt  y  coordinates  are  the  coordinates  of 
the  shadow  axis"  intersection  with  the  fundamental  plane. 


C.3  PENUMBRAL  AND  UMBRAL  RADII  ON  1  HE  FUNDAM ENTAL P LANE 

In  Figure  (.  ?  two  shadow  cones  have  been  drawn  in  cross  section,  the  interior 
tangent  cone  defining  the  penumbra  :  nd  having  the  vertex  Vp  a  distance  c,  above 
the  fundamental  plane  and  the  exterior  tangein  eor.e  defining  the  umbra  and 
having  a  vertex  V  a  distance  c2  above  the  fundamental  plane.  We  may  write 
immediately  that 


sin  /, 

00  T  ci^ 

G 

(C.6) 

sin  f2 

_  d©  - -dfc 

G 

(C.7) 

where  /  and  f7  are  the  penumbral  and  umbral  cone  half-angles,  respectively,  and 
where  dG  and  dv  are  the  linear  semi-diameters  of  the  sun  and  moon,  respectively. 
Rut,  we  not  ’  from  (C.2)  and  (A. 22)  th3t 


G  =  grG  -  g 


-43—  =  ,  _i?_  R, 
sin  7f0  sin  rr0 


(C.8) 


.0111  the  definition  of  the  solar  semi-diameter  that  sin  Sq=  dG/rQ,  and  fur- 
ibc  from  (A. 23)  that 


d@  =r@  sir.  S&.  -■  (sin  S0/sin  v0)  ae;  (C.9) 

hence,  (C.6)  and  (C.7)  can  be  written  in  the  usual  form  as 

sin/,  =  (sin  S0  +  k  sin  TT(,)/gR  (C.10) 

sin/2  =  (sin  S0  -ksin7r0)/gK  (C.  11) 

where  k  =  dt/affi  is  an  adopted  constant. 


From  AVpBC  we  may  write  that 


c,  =  z  +  dt  cosec  /, 
and  similarly  from  AVuB'C' 

<h  =  z  -  d4  cosec  f2. 


hence. 


and 


ct  =  (z  +  k  cosec  /  >  a€ 


Cj  =  (z  -  k  cosec /j); 
where,  from  (C.4)*  and  (C.5)  respectively. 


z  =  za^ 


and 


(C.12) 


(C.13) 

(C.14) 

(Cl  5) 


2  =  Jsln  sin  d  -  cos  6V  cos  d  cos  (o^  -  a)]  /sin  jtv. 

Thus,  the  penumbral  radius  on  the  fundamental  plane  can  be  found  from 

S,  =  c,  tan/,  (C.16) 

and  the  umbral  radius  can  be  found  from 

6j  =  oj  tan  /j .  (C.17) 

The  sign  convention  introduced  in  (C.  15)  is  such  that  £3  is  negative  for  total 
eclipses  and  positive  for  annular  eclipses. 

C.4  SUMMARY  -  THE  BESSELI AN  ELEMENTS 

In  the  developments  of  the  prior  sections  the  geometric  position  of  the 
shadow  axis  and  radii  of  the  penumbral  and  umbral  cones  on  the  fundamental 
plane  have  been  described.  The  elements  appropriate  to  this  description  are  the 
Besselian  elements,  x,  y,  sin  d,  cos  d,  £, ,  £3  and  n,  the  ephemeris  hour  angle  of  the 
cJhadow  axis.  This  hour  angle  is  used  in  place  of  the  a,  the  right  ascension  of  the 
shadow  axis,  and  is  calculated  from 

M  =  E.S.T.-a  (C.18) 

where  the  ephemeris  sidereal  time  E.S.T.  is  given  by  (B,3). 

*  We  note  thet  a4.  the  earth's  equatorial  radius  is  generally  set  equal  to  unity.  In  such  instances, 
of  course,  (C.4)  immediately  reduces  to  the  trivial  identity  z  -  z. 


These  elements  and  their  hourly  variations,  which  are  independent  of  posi¬ 
tion  on  the  spheroid,  are  tabulated  in  the  major  almanacs  and  the  major  eclipse 
canons,  such  as  (4) 

C.5  COORDINATES  OF  THE  OBSERVER  RELATIVE  TO 
THE  FUNDAMENTAL  PLANE 

We  will  consider  a  spheroid  of  flattening  f  (or,  equivalently,  of  ellipticity  or 
first  eccentricity  e>*  centered  at  the  geocenter  and  having  an  equatorial  radius  as 
where 

<is  —  3f  a^  ( C .  1 9 ) 

with  3f  an  arbitrary  constant;  positions  will  be  specified  on  the  spheroidal  surface 
by  assignment  of  longitude  (West  taken  positive)  and  geodetic  latitude.  Hence,  by 
employing  the  same  transformation  used  to  effect  (C.4)  and  (C.5),  an  observer 
at  ephemeris  longitude  and  latitude  0  will  have  the  fundamental  plane 
coordinates 

£  *  Ip  cos  <p'  sin  (h)  ]  K  a® 

tj  =  [p  sin  0'cos  d  -  p  cos  <p‘  sin  d  cos  (h)  ]  Jfa® 

J  =  (p  sin  0'  sin  d  +  p  cos  0'  cos  d  cos  (h)  ]  Jfa®  (C.20) 

The  hour  angle  (h)  is  given  by 

0=P-Ae  (C.21 ) 

where  the  ephemeris  longitude  Ae  is  related  to  the  longitude  A  by  means  of  the 
relation 


A,  =  A  +  1.002738  AT.  (C.22) 

AT  is  the  difference  (either  extrapolated,  in  the  case  of  the  future  or  distant  past, 
or  measured,  in  the  case  of  the  immediate  past)  between  Ephemeris  Time  and 
Universal  Time.  The  geocentric  latitude  0'  and  geocentric  distance  of  the  observer 
p  are  related  to  the  geodetic  latitude  0  through 

p  sin  0'  =  S  sin  0;  P  cos  0'  =  C  cos  <t>  (C.23) 


*  We  note  from  the  definitions  /  =  (a  —  b)/a  and  e  =  [1  —  b2  /  a2] !j,  where  a  is  the  equatorial 
radius  end  b  the  polar  radius  of  the  spheroid,  that  e3  =  2f  -  f3 . 
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where 


S  =  (l  -f)2  C  C=  [cos2  0  +  (l  -/)2  sin2  0]  ~w  (C.24) 

or,  equivalently 

S  =  ( 1  -  e2  )C;  C  =  ( 1  -  e2  sin2  0)  ~  w  (c. 25) 

C..6  PENUMBRAL  AND  UMBRAL  RADII  ON  THE  OBSERVER'S  PLANE 

We  now  construct  a  plane  through  the  observer  parallel  toandat  a  distance  f 
above  the  fundamental  plane.  Referring  once  again  to  Figure  C-2,  it  is  seen  that 
the  penumbral  and  umbral  radii  on  this  plane  are  given  by 


L,  =(c,  -  f)  tan  /(  =  8,  -  f  tan /, 


(C.26) 


and 


L2  =  (c2  -  f)  tan  /j  =  82  -  f  tan  /2 


(C.27) 


respectively. 

With  these  two  equations  and  some  additional  analysis,  it  becomes  possible 
to  calculate  all  circumstances  of  the  eclipse  appropriate  to  the  locale  of  a  specified 
observer;  these  calculations  will  be  discussed  extensively  in  a  subsequent  appendix. 
But,  if  one  wishes  to  predict  the  observation  site  or  sites  appropnate  for  a  given 
circumstance  -  such  as  the  outline  of  the  umbra  on  the  spheroid  -  then,  the 
spheroidal  flattening  introduces  the  complication  that,  in  order  tc  calculate  L2 
from  (C.27),  for  example,  f  must  be  known  which  requires,  in  turn,  that  0'  and  0, 
a  desired  end  result,  must  be  known.  This  difficulty  is  minimized  by  the  following 
treatment  of  Bessel. 

C.7  INTRODUCTION  OF  THE  AUXILIARY  ELEMENTS 

The  development  is  initiated  by  introduction  of  the  parametric  latitude  0, , 
through  the  relation 

cos  0t  =  P  cos  0’  =  cos  0/(1-  e2  sin2  0)^  (C.28) 

from  which,  using  (C.25), 

sin  0,  =  [1  -  cos2  0,  !!/’  =  (1  -  e2 )"  w  Psin  0'  .  (C.29) 
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Consequently ,  (C  20)  becomes 


£  =  |cos0,  sin  ©  ] 

tj  =  ((I  —  c3  )Vi  sin  0,  cos  d  -  cos  0,  sin  U  cos  ©  J  -Ka^ 
f  =  ((1  -  e2  )w  sin  0(  sin  d  +  cos  0,  cos  d  cos  ©  J  Ufa®  (C.30) 

We  now  make  the  following  substitutions  in  r(, 

sin  d  -P,  sin  d,  ;(1  -  e2  y1  cos  d  =  P,  cos  d,  (C.31 ) 

and  the  following  substitutions  in  f, 

cos  d  =  p2  cos  dj ;  (1  -  e2  y’  sin  d  =  P3  sin  d2  (0.32) 

so  that  (C.30)  now  becomes 

(  =  t cos 0]  sin  ©  ]3Cae 

t?,=  VIP  i  =[sin0,  cosd,  -  cos  0,  sin  d,  cos©  1  JCag, 

f  =  Pj  l  sin  0,  sindj  +cos0,  cos  d3  cos  ©  lKa®.  (C.33) 

Next,  we  define  the  variable  J,  from 

f,  s  (Xa®)2  -*2  -»)2  (C.34) 

a  useful  form  of  which  -  after  substitution  of  £  and  i?  from  (C.33)  and  some 
reduction  -  is 

f,  -  (sin  0,  sin  d,  +  cos  0,  cosd,  cos  (ri)  1  Jta^,  (C.35) 

Further,  we  note  from  (C.33)  that  f,  and  rj|  may  be  readily  manipulated  to  yield, 


3Ca®(sin0,)  =  r?,  cosd,  +  J,  sind, 

and 

JCa®  (  cos  0,  cos  ©  )  =  -Tj,  sin  d,  +  f ,  cos  d, 
which  upon  insertion  into  f  of  (C.33)  yields 


(C.36) 

(C.37) 


?  =  Pi  lfiCOs(d,  -  dj)  -  j?,  sin  (d,  -  d3)).  (C.38) 
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This  torm  ot  f  lias  the  useful  property  that  the  assignment  of  specific  values 
to  £  and  tj  -  hoth  of  which  arc  in  the  fundamental  plane- allows  the  immediate 
calculation  of  f ,  by  (C.35)  and  thus  the  third  coordinate  f  by  (C.38). 


f  inally,  from  the  first  equation  of  (C.33)  writ*,'-,  in  the  form 
Tf a,^  (cos  <*>,  sin  (h)  )  =  £ 
and  from  (t'.3b)  and  (\37),  it  follows  that 

t 


tan  (h)  = 


f,  cos  d,  -  r},  sin  d, 


and 


sin  0,  =  — L  (J,  sin  d|  +  rjt  cos  d, ). 

(Ifa^) 


(0.3d) 


(C.40) 


(0.41) 


Thus,  if  definite  values  of  £  and  tj  are  available  from  the  application  of  specific 
geometric  conditions,  {.  tj,  and  ?,  can  be  calculated  and.  hence  (h)  land  thus  X) 
and  <<>|  (and  thus  <t> )  can  be  obtained  from  (('.40)  and  (C.41 ),  respectively 


These  calculations,  and  the  calculation  for  f  in  (C  .38),  fully  incorporate  the 
flattening  and  are  based  upon  the  auxiliary  elements  P, .  Pj,  sin  a,,  sin  d3. 
sin  (d,  -  dj)  and  cos  (d,  -  d2 )  which,  after  manipulating  (C.31)  and  (C.32), 
follow  directly  from  the  Besselian  element  d  and  the  ellipticity  of  the  spheroid  e  . 
in  the  form 

P,  =  (1  -  e1  cos2  d)h 
Pi  =  ( 1  -  e}  sin2  d)'' 


sin  d, 

=  sin  d/P , 

cos  d, 

=  (1  -  e2  )1'  cos  d/P, 

sin  (d, 

-d,) 

=  e2  sin  d  cos  d/P,  P 

cos  (d 

,  ~d,) 

=  0  -  e2  )1’ / P,  P,  . 

(C.42) 
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C.8  OUTLINE  OF  THE  SHADOW 


In  general,  eclipse  phenomena  are  described  by  specifying  tire  dis'ance  and 
position  angle  0  (measured  cast  from  north)  of  the  shadow  axis  with  respect  to 
the  “observation''  point  in  question.  Tims,  in  the  plane  drawn  through  the 
“observer"  and  maintained  perpendicular  to  the  shadow  axis,  we  can  write  thai 
the  mi- plane  coni|K>nents  of  the  observer  coordinates  (£.  rj)  arc  related,  te>  the 
shadow  axis  coordinates  t\,  y)  In  means  of 

i  =  x  m  sin  Q 
and 

t /  -  y  m  cos  Q  (0.44) 

where 

m1  =  (x  <)*’  +  (y  tj);  ((.  ,-iS) 

Consequently,  it  follows  immediately  that  the  outline  of  the  umbra*  on  the 
spheroidal  surface  at  a  specific  time  is  given  by  those  (£.tj)  pairs  which  circum¬ 
scribe  (x.y)  at  the  distance  m  =  L2 . 

However,  before  allowing  Q  to  assume  alt  values  from  0°  to  M<0'  in  calcu¬ 
lating  the  sliadow  outline,  it  is  useful  to  determine  whether  the  entire  umbra!  cone 
intersects  the  spheroidal  surface  or,  equivalently,  whether  the  entire  utnbral  circle 
falls  within  a  circle  of  (nearly  constant)  radius  K'a.j,  on  the  fundamental  plane. 

Figure  C-d  portrays  the  situation  when  a  portion  of  the  umbra  1  circle  falls 
outside  this  circle. 

If  we  now  set  x  =  7  sin  in  and  y  =  7  cos  m  (where  M  is  the  position  angle  of 
the  shadow  axis)  and  note  that  at  the  extreme  points  Lj  =  Vj  (since  these  points 
are  in  the  fundamental  plane  and  f  =  0),  then 

t.NXr,)-'  =  T  :  +  V’j  :  y  V':  cos  (Q  -  M) 

or  ... 

T*  -t-c;  -(Nar¬ 
cos  (0  M)  = -  «.  .4(0 

;v;  7 


The  outline  ot  the  penumbra  can  be  calculated  by  replacing  V ;  and  L;  by  V,  and  l 
respectively 


where  7 3  -  x3  +  y3 .  When  lhe  evaluation  of  the  right-hand  member  of  ((  40) 
leads  to  the  condition  that  cos  (Q  -  M)  <  1,  the  two  consequent  values  of  Q 
define  two  circumferential  segments:  one  within  the  fundamental  plane  circle  and 
(running)  values  of  Q  producing  J  >  0,  the  other  outside  the  fundamental  plane 
circle  and  (running)  values  of  Q  leading  to  f  <0.  When  (t'.4(>)  results  in 
cos  (Q  -  M )  >  1 ,  the  entire  shadow  lies  within  the  fundamental  plane  circle  and 
all  values  of  Q  from  0°  to  360°  are  peimissiblc. 

Assuming  that  the  appropriate  range  of  values  lias  been  established,  t l*c 
shadow  outline  can  now  be  calculated  by  tin-  following  iterative  procedure  for  each 
of  the  permissible  Q. 


Step  1:  Assume  La  =  li}  and  calculate 

£,('»  =  x  —  Ej  sinQ 
’t\0)  =  (y  -  Ej  cosQ)/p, 

$10>  *  [(Jfa®)1  -  ({<n,):  -  (rjVH)J  lv*  (C  -47) 

Step  11:  From  Vi01  and  ft01  calculate  f<0)  from 

?0)  =  p2  |fl0>  cos  (d,  -  d2 )  -  ij\°‘  sin  (d,  -  da  )1  (C.48) 

and  LEui  from 

U°»  =  fij  -  r("  tan  /,  (C'  49) 

Step  111:  Using  this  value  of  La  calculate 

£(1)  =  x  -  L^0>  sin  Q 
rjU)  =  (y  -  LEn)  cos  Q )/p, 

r,,l)  =  l(Ha,j.)J  -({'"l2  -(T?l"):  )V;  ((.50) 

Step  IV:  Assuming  that  an  additional  iteration  of  Step  11  is 

not  required  by  the  results  of  (C.50),  use  (('.50)  to 
calculate  -  following  (C  40)  and  (C.41 )  - 


(an  (u) 


f1<1>  cos d,  -  r?!1 1  sin  d, 


(('.51) 
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and 


sin  <P, 


Si 


3fa, 


(C  52) 


l9> 


which  leads,  following  (C.21 )  to 

\  =  »  -  ®  (C-53) 

and,  following  (C  23),  (C.25S,  (C.28)  and  (C.29)  to 

tantf  =  <1  -c1)'*  tan  0,  (C.54) 


C.9  THECENTRAL  LINE 

Sirce  the  central  line  is  the  locus  of  points  of  intersection  of  the  shadow  axis 
with  the  spheroid,  it  follows  immediately  that  the  results  of  the  previous  section 
may  be  utilized  by  shrinking  the  shadow  radius  to  zero.  Thus,  the  equations  of 
(C.50)  become 


{  =  x 

t?i  =  y/p,  =  yi 

f,  =  (CKa*)5  -xJ  -y;)*  (C.55) 

which  leads  immediately,  upon  substitution  in  (C. 51)  and  (C.52)  to 


tan  (h )  -  - 

v  j  ..  J 

Si  Uj  —  V  j  MM  Uj 

ana 

sin  d,  +  y,  cos  d, 
sin  0,  = - 

where  again  from  (C.54) 

tan  <t>  =  (1  -  eJ )'  w  tan  0,  . 


(C.56) 


(C.57) 


(C.58) 
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If  (('  4)  and  ((’.34)  ere  employed,  then  (C  .50)  and  ((’.57)  aie  transformed  to 


tan  (m)=- 


r  * 2  +  v;  1 

Vj 

K 

1 

cos  d,  -  >  ,  sin  d 

(f.59) 


and 


i_ 

j 

1 

j 

a3  +  y? 

sin  0,  =  — 

l  X 

1  1 

1  X 

{ 

L  x2  } 

sind,  +  y,  cosd,  I  ((’.60) 


which  is  the  form  employed  by  Lewis*  in  calculating  the  ionospheric  central  line. 

C.10  NORTHERN  AND  SOUTHERN  LIMITS  OF  THE  SHADOW  - 
ADDITIONAL  AUXILIARY  ELEMENTS 

From  the  discussion  immediately  following  ((7.45),  it  is  readily  seen  that  the 
examination  of  the  time  variation  of  the  observer/shadow-cone  separation  mJ  -  L3 
(wheye  L  is  either  the  pcnumbral  or  umbral  radius)  can  be  employed  either  to 
calculate  circumstances  given  the  observation  location  or,  conversely,  to  calculate 
the  observer  location  given  the  circumstances.  Here  wc  are  concerned  with  the 
latter  type  and,  more  specifically,  with  calculating  the  northern  and  southern 
limits  of  the  shadow  directly  rather  than  by  inferring  these  values  from  the  com 
plete  shadow  outline  calculated  in  Section  C.8. 

In  Figure  C-4,  a  time  variation  of  nr3  -  L1  has  been  hypothesized  for  two 
locations.  At  location  A  the  particular  phase  of  the  eclipse  under  consideration 
(partial  for  L  -*  L, ,  total  for  L  -+  L2)  begins  at  t(,)  and  ends  at  t(f\  i,e.,  at  the 
two  solutions  ot  ni3  -  L3  =  0.  At  location  B,  on  the  other  hand,  the  eclipse  both 
begins  and  ends  at  the  same  time,  an  occurrence  identified  by  the  coalescence  of 
the  two  roots  of  m3  -  L3  =  0  to  tR  and  the  presence  of  the  horizontal  tangent  at 
tB  Hence,  at  tB  we  liavc  that 

(x  -  $)3  +  (y  -  r?)3  -  (8  -  f  tan  /)3  =  0  (C.61 ) 

and,  using  dots  to  signify  diflercntiation  with  respect  to  time,  that 

(x  -  {)  (x  -  £)  +  (y  -  rj)  (y  -  r?)  -  (5  -  f  tan  /)((-{  tan  f)  =  0  (C  .62) 
where  /,  given  by  (0  i  0)  or  (C.  i  i ),  is  sensibly  constant. 


Lewis,  op.  cit. 
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FIGURE  C— 4  HYPOTHETICAL  TIME  VARIATION  OF  m5  -  L1  FOR  TWO  LOCATIONS 

We  now  note  that  differentiation  of  (C.20)  leads,  after  a  modest  reduction,  to 

£  =  n  (-  t)  sin  d  +  f  cos  d) 

•  •  • 
tj  =  p£sind-df 

f  =  p£  cosd+drj  (C.63) 

(j.i  and  d  are  presented  as  hourly  variations  in  the  literature)  which  become,  upon 
substitution  of  (C.43  and  (C.44) 

£  =  P  (-  y  sin  d  +  f  cos  d  +  in  sin  d  cos  Q) 

V  =  p  (x  sin  d  -  m  sin  d  Q)  -  d  f 

f  -  p(  —  x'osd  +  m  cos  d  sin  Q)  +  d  (y  -  m  cos  Q).  'C.64) 

Hence,  (C.62;  becomes,  upon  substitution  of  (C.64),  and  upon  noting  that  at 
t  =  t„,  m  =  L  =  8  -  f  tan  f. 
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sin  Q  (x  +4  y  sin  d  +pt  2  tau  /  cos  d  -  £  cos  d  (1  +  tan2  /)]  + 

+  cos  Q  (y  -  p  x  sin  d  +  d  f  -  d  tan  /  (8  -  f  tan  f) )  + 

+  [  8  —  p  x  cos  d  tan  /  +  d  tan  /  y]  =  0.  (C.65) 

or,  upon  defining  the  auxiliary  elements 
a  =  -S-/ix  tan  /  cos  d 
b  =  y  +  fi  x  sin  d 

c  =  x  +/i  y  sin  d  +/i  8  tan  f  cos  d  (C.66) 

(where  a,  and  c,  are  penumbral  auxiliaries;  a2  and  c2  are  uir.bral  auxiliaries) 
sin  Q  (c  —  ^  f  cos  d  ( i  +  tan3  f) ]  + 

eosQ  (- b  +d  J  -  d  tan  /  (8  -  £  tan  /)]  +a  +d  tan  i  y  =  0.  (C.67) 

However,  following  Cha  ivenet,  terms  involving  tan2  J  and  d  tan  /  will  be  neglected 

so  that  we  have 

tanQ=  irliLz»J5£0_ .  (C.68) 

c  -  M  t  os  d 

as  the  condition  for  the  northern  and  southern  limits.  For  umbral  limits,  cos  Q  is 
positive  for  the  northern  limit  of  a  total  eclipse  and  the  southern  limit  of  an  annu¬ 
lar  eclipse;  cos  Q  is  negative  for  the  southern  limit  of  a  total  eclipse  and  the 
northern  limit  of  an  annular  eclipse.  F  .r  penumbral  limits,  cos  Q  is  positive  for 
the  southern  limit  and  negative  for  the  northern  limit. 

The  calculations  for  the  limits  are  similar  to  those  of  the  sliadow  outline 

Step  I:  Assume  L?  =  £j  (or  Li  =  8, )  and  calculate 

tan  Qt0)  =  (C.69) 

c 


from  which 


t<°>  =  x-£,  sin  Q!0) 

=  :y-C2  cos  Q(0,)/Pi 

fjO)  =  -(i,0))3  -(rj,<0))J]w  (C.70) 

Step  II:  From  »}[0)  and  fj01,  calculate  f<0) 

f10’  -Pi  Ui0)  cos(d,  -  dj)-^0)  sin(d,  -d2)J  (C.71) 

and  thence 

b  -  u  ^0)  -  a^  sec  Q(0) 
tan  Q(  u  = - - 

Cj  —  M  f <0>  cos  d  (C.72) 

and 

Un  =  fij  -  ^0)  tan/j  (C.73) 

Step  III:  Using  these  values  of  QU)  and  Ljl),  calculate 

£<l>  =  x  -  Lj(1>  sin  Q(1) 
tj(>)  =(y_L^i)  cosQ(,))/Pi 

=(3faffl)3  -«<n y  -(„|»>)*]*  (C.74) 

Step  IV :  Assuming  that  an  additional  iteration  of  Step  1!  is  net 
required  by  the  results  of  fC.74),  proceed  as  in  (C.51 ), 
(C.52),  and  (C.53). 


C.11  DURATION  ON  CENTRAL  LINE 

Given  an  observation  point  £0,  t?0  and  a  center  line  point  x0.  y0  both  evalu¬ 
ated  at  the  time  T0,  the  two  components  of  the  position  of  the  shadow  axis 

fplntivp  to  t hie  nhcprvatirtn  nA«nt  Ka  AvnA^^r.4  r:- — * - 1 —  _*  4l  .  . 

- -  _  vVw.  .ui.vi.  KVi.iv  mu  i/v  vApomiCu  »o  tudi  uiuci  ai  me  tunc 

T0  +  r  as 


(x  -  i)  =  [x0  +  x  r  ...)  -  [fu  +  £  r  +  ...)  =  m  sin  Q 


(C.75) 
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I 


I 

| 


(y  -  »j)  =  (y0  ry  t  +  ...)  -  li)0  +tj  r  =  m  cosQ 


(C.76) 


or,  upon  squaring  and  adding  the  components  to  yield  the  relative  separation, 
m3  =  (x0  -  i0)J  +  (yo  -  no)2  +  2[(x0  -  £0)  (x  -  t)  +  ty0  -  no)  <y  -  n))*  + 


+  l(x  -  +  (y  -  n)2 1  +  ... 


(C.77) 


We  will  now  place  the  observer  on  the  center  line,  such  that  £0  -  x0  and 
n o  -  y0  and  such  that  the  time  r  =  T0  -  Ti  will  be  taken  as  the  beginning  of  the 
umbral  eclipse  when  m  -  L2 . 

Thus  (C.77)  becomes 


Lj  =  l(x-£)3  +(y  -n)2)  (To  -T,)3 


(C.78) 


and  hence,  the  semi-duration  s  =  T0  —  T;  of  the  eclipse  is  immediately  given  by 
the  solution  of  (C.80)  in  the  form 


(C.79) 


where  n,  the  speed  of  the  shadow  relative  to  the  observer,  is  given  by 


n  =  [(x  -  £)2  +  (y  -  n)2  1  . 


(C.80) 


Lj  is  assumed  to  be  constant  throughout  the  calculation  and  is  evaluated  from 


U  =  -  5"o  tan  {, 


(C.81) 


in  which 


fo  -  Pi  (ii  cos  (d,  -  d2)  -  ■»),  sin  (d:  -  d2)l 


(C.82) 


where  f,  and  tj,  are  evaluated  from  (C.55)  with  x0  and  y0  replacing  x  and  y. 
T)  can  be  evaluated  by  noting  that  since  £  =  x  and  t)  =  y,  (C.63)  can  be  manipulated 
to  yield 


x-£=x  +  Mysind-pf  cos  d 


(C.83) 


y-’J-y-^xsind+df 


(C.84) 


or,  by  employing  the  auxiliary  elements  of  (C.66), 


(C.85) 


x  —  £  =  c  -  m  cos  d  (f  +  C  tan  f) 
y— ij=-b+di 


(C.86) 


which,  consistent  with  the  level  of  approximation  of  (C.75)  through  (C.80), 
reduce  to  the  approximate  relations, 


c  -  n  f  cos  d 


(C.87) 


(C.88) 


APPENDIX  D 


LOCAL  CIRCUMSTANCES 


DO  INTRODUCTION 

As  pointed  out  in  Section  C.5  of  the  previous  appendix,  it  is  possible  to  pre¬ 
dict  all  circumstances  of  an  eclipse  at  a  locale  once  the  coordinates  of  that  locale 
have  been  specified.  The  purpose  of  this  appendix  is  to  develop  the  basis  of  algo¬ 
rithms  employed  to  calculate  these  circumstances,  to  develop  differentia!  correc¬ 
tion  piocedures  for  adjusting  these  circumstances  to  modest  departures  from  the 
specified  locale  and  to  develop  the  calculations  of  other  topocentric  parameters  of 
potential  interest  to  solar,*  atmospheric  and  ionospheric  research. 

D.1  PRELIMINARY  CONSIDERATIONS 

Following  (C.20)  we  may  write  immediately  for  a  site  with  geodetic  latitude 
0,  longitude  X  (West  taken  positive)  and  height  h  (in  meters)  above  the  spheroid 
that 


£  =  p  cos  0'  sin  (h 

t)  =  p  sin  </>'  cos  d  -  p  cos  0'  sin  d  cos  (h) 

f  =  p  sin  <t>"  sin  d  +  p  cos  0'  cos  d  cos  0  ( D.  1 ) 

in  which  the  equatorial  radius  of  the  earth  a€,  will  be  assumed  henceforth  as  the 
unit  of  distance  and  where 


p  sin  0'  =  (S  +  H)  sin  0;  p  cos  0'  =  <C  +  H)  cos  0  (D.2) 

with  H  =  0.15678503  x  10'6h;**  S  and  C  are  given  by  (C.24)  and  (C.25).  The 
hour  angle  (0  is,  following  (C.21 )  and  (C.22), 

(h)  =  p  -X  -  1.002738  AT  (D.3) 


’Appendix  E  is  devoted  to  the  development  of  various  coordinate  systems  (and  their  inter¬ 
relationships)  for  use  ;n  observation  of  optical  and  radio  phenomena  on  and  above  the  solar 
surface. 

**We  recognize,  witn  reference  to  (C.20),  the  difference  between  the  change  of  scale  implied 
by  Jfa$  of  (C.20), on  the  one  hand,  and  the  addition  of  the  height  term  in  D.2  normal  to  the 
spheroid,  if  added  at  each  point  of  the  spheroid,  on  the  other.  However,  for  a  scale  change 
(or  height  addition)  of  O.la^the  difference  between  the  two  approaches  leads  to  a  difference 
of  2km  at  the  poles  and  will  be  neglprted. 
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where  n,  the  ephemcris  hour  angle  of  the  shadow  axis,  is  given  by  (C.  1 8)  and  AT, 
is  the  difference  of  Ephemeris  Time  and  Universal  Time  for  the  year  in  question. 

If  these  coordinates  and  their  time  variations,  given  by  (C.63)  as 


£  =  v  ( -  r?  sin  d  +  f  cos  d) 
j?=%i£  sind-df 
?  »  £  cosd  +d  tj, 

are  combined  as 

u  =  x  -  £;  v  -  y  -  rj;  m2  =  u2  +  v2 


and 


u  =x-£;v=y-T?;n2=u2+v2  , 


(D.4) 


(D.5) 


(D  6) 


then  (C.77),  the  first  order  expression  relating  the  distance  m  between  the  shadow 
axis  (x,y)  and  an  observation  site  (£,  f?)  at  the  time  T0  +  r  to  the  shadow  axisand 

observation  site  positions  (x<) ,  y0 ;  £0 ,  »?<>)  and  velocities  (x,  y;  £,  fj)  at  the  time  T0, 
becomes 


m2  =  m2  +  2(uo  u  +  v0  v)  r  +  n2  r2  . 
This  equation  has  the  solution 

r  =  _  _D  +  _D_ 


»-  ~  (mo 


(D.7) 


(D.8) 


where  D  =  u*,  u  +  v0  v,  or,  following  rearrangement  of  the  second  term 


,  __L_ 

(  u0  v  -  v0  u 

VI 

m* 

l 

)  J 

r  =  _  _D  ±  JH 

n2  n 


which  leads  to  the  to.  m  commonly  encountered  in  the  literature 
D 


T  =  -  n 2 


m  _ , 

-  U/5  V' 

n 


(D.9) 


where  sin  =  (u0  v  -  v0uVn  m. 
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Thus,  to  determine  the  time  at  which  ni  assumes  a  specified  value,  one  esti¬ 
mates  an  Ephemeris  Time  T0  “reasonably"  close  to  the  desired  time,  evaluates  all 
component  quantities  of  (D.5)  and  (D.6)  and  solves  -  after  insertion  of  the  speci¬ 
fied  value  of  m  -  either  (I). 8)  or  (D.9)  for  the  correction  r  to  the  estimated  time. 

If  great  accuracy  is  required,  a  second  iteration  may  be  required  in  which 
(D.5)  and  (D.6)  are  re-evaluated  at  T0  +  r  -  T^,  and  the  correction  r'  to  T|,  is  cal¬ 
culated  from  (D.8)  or  (D.9). 


Similarly,  by  differentiating  (D. 7)  and  rearranging  the  result  to  yield 


(D.  10) 


one  can  determine  the  increment  t  to  the  time  T0  at  which  m  assumes  a  specified 
value,  all  other  “component”  quantities  of  (D.  10)  being  evaluated  at  the  “reason¬ 
ably”  close  time  T0 . 

D.2  TIME  OF  GREATEST  PHASE 


The  time  of  greatest  phase  occurs  when  the  eclipse  magnitude  M,  =  ( L,  —  m)' 
(L,  +  L2),  derived  in  Section  D.7,  is  precisely  a  maximum.  But  since  the  shadow 
radii  variations  L,  =*  L2  =  L  are  small,*  the  time  of  greatest  phase  is  usually  taken 
as  the  time  at  which  the  distance  m  between  the  shadow  axis  and  the  observer  is 
minimum,  or  when  m  -  0. 


If  this  specified  value  of  m  is  substituted  in  (D.10),  then  the  correction  t  to 
the  "reasonable”  time  of  greatest  phase  T0  is  immediately  given  by 


T 


(D.l  I) 


which  leads  to  the  interpretation  of  the  first  terms  in  both  (D.8)  and  (D  9)  as  that 
correction  to  the  assumed  time  T0  which  effects  the  minimum  separation  between 
the  shadow  axis  and  the  observer. 


If,  on  the  other  hand,  accuracies  of  tens  of  milliseconds  are  required,  then 
the  maximum  magnitude  definition  employed  by  Gossner** 


M. 


d  /  L.  —  m  \ 

if  lTTr:ru 


(D.l  2) 


*  cf.  (D.39)  and  associated  fooinote. 

**  Gossner,  S.  D.,  "A  Correction  to  the  Time  ol  Maximum  Obscuration  in  Solar  Eclipses,' 
A,  J.  60,  383  U955). 
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which  can  be  written,  following  substitution  of  L,  =  L2  =  L,  in  the  form, 


m 


LLit  *  i-r)L)  C  L  i  -  m)  (L|  +  L3 )  ^  (La  —  L|  +  2m) 


(L,  +  L-) 


Li  +  Li 


(D.  1 3) 


leads,  upon  substitution  of  (D.13)  in  (D.10)  to  Gossner's  Eq  (5) 


/  l  -i  L,  +  2m 
\  L,  +  La 


(D.14) 


which  reduces  immediately  to  (D.  1 1 )  when  shadow  radii  variations  are  neglected. 


D  3  PENUMBRAL  (FIRST  AND  FOU RTH)  CONTACTS  AND  DURATION 


In  the  case  of  the  penumbral  (first  and  fourth)  contacts  we  set  m  =  L,  in 
(D.9),  following  the  choice  of  some  initial  time  T0,  and  write  immediately  that 


in  which 


Up  V  -  Vp  u 
n 


(D.15) 


(D.16) 


the  negative  sign  is  selected  in  calculating  the  correction  to  the  Ephemeris  Time  T0 
for  first  contact  (immersion)  and  the  positive  sign  is  selected  for  fourth  contact 
(emersion).  Since  L,  loos  J/p  l/n  is  either  added  to  or  subtracted  from  the  time  of 
maximum  eclipse  T0  -  D/nJ  to  give  either  the  first  or  fourth  contact,  then  assum¬ 
ing  T0  has  been  selected  to  be  reasonably  close  to  maximum,  Lj  loos  l/n  can  fee 
regarded  as  the  semi-duration  of  the  partiai  phase. 


Because  the  times  of  penumbral  contacts  are  generally  not  required  with 
great  accuracy,  one  can  ,-and  the  square  root  of  (D.8)  and  keep  only  the  first 
term,  thus  yielding  as  the  approximate  correction 


r 


p 


±  <L?  ~nio) 

2D 


(D.17) 


where  again  the  negative  and  positive  signs  are  selected  for  first  and  fourth  con¬ 
tacts,  respectively. 


D.4  UMBRAL  (SECOND  AND  THIRD)  CONTACTS  AND  DURATION 


l  or  the  second  and  third  (umbral)  contacts  m  =  Lj  and  (D.9)  becomes 


D  4  1,  . 

’■  ■  - 1? 1  -r  1  c“  * 


in  which 


sin  y 


n„  v  -  v„  11 


(D.18) 


(D.19) 


The  positive  sign  is  selected  either  lor  second  contact  (immersion)  in  the  case 
of  a  total  eclipse  (L3  <■  0)  or  third  contact  (emersion)  for  an  annular  eclipse,  the 
negative  sign  is  selected  either  tor  third  contact  (emersion)  for  a  total  eclipse  or 
second  contact  (immersion)  for  an  annular  eclipse  Because  of  accuracy  require¬ 
ments  a  second  iteration  is  often  required. 


The  semi-duration  of  the  umbral  please  is  given  immediately  as 
l-i 


Su  »  n 


COS  l)/ 


(D.:o) 


which  reduces,  as  expected,  to  the  center  line  value  of  l.2/n  given  by  (C.81 )  when 
Uo  =  v0  =  0  and  cos  t!/u  =  1 . 


D.5  POSITION  ANGLES 


The  angle  Q  which  can  be  calculated  —  following  <1. 43),  tC.44),  (D.5)  and 
(D.6)  -  from 


tan  Q  =  y  (D.2I) 

is  at  all  times  the  topocentric  position  angle  of  the  moon’s  center  measured  east¬ 
ward  from  the  north  point  of  the  solar  limb.*  When  u  and  v  aie  evaluated  at  the 
contact  times  of  (D.15)  and  (D.  18),  Q  is  the  angle  of  the  points  of  contact  of 
the  apparent  disc  of  the  moon  with  the  apparent  disc  of  the  sun. 

In  instances  when  altitude-azimuth  instrumentation  is  being  employed,  the 
position  angle  with  respect  to  the  vertex  (defined  as  that  point  of  the  solar  limb 
nearest  the  zenith)  may  prove  more  useful.  This  angle  is  determined  by  noting  that 


*  The  north  point  of  the  sun  is  readily  determined  observationally  as  that  point  on  the  apparent 
solar  disc  tangent  to  the  disc's  apparent  (diurnal)  motion  since  it  is  the  intersection  point  of 
the  northern  limb  and  the  hour  circle  through  the  solar  center. 
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the  distance  from  the  nortli  point  of  the  sol.,,  limb  to  the  vertex  is  equal  to  the 
parallactic  angle  q  which  can  be  gotten  from  the  equation 


tan  q  - 


cos  0  sin  0 


sin  Q  cos  d  -  cos  <p  sin  d  cos  0 


(0.22) 


or,  since  cos  $  =%  p  cos  <?'  and  sin  <p  ^  P  sin  in  (D.l),  fiom  the  approximate 

relation 


,  _  t 

tan  q  ^  — 
9 


(D.23) 


Hence,  the  position  angle  V,  as  measured  from  the  vertex,  is  given  by 

V  =  Q  -  q  .  (0.24) 


D.C  DIFFERENTIAL  CORRECTIONS  TO  LOCAL  CIRCUMSTANCES 

In  this  section  all  differential  corrections  requires  to  adjust  calculated  local 
circumstances  to  modest  cltanges  in  locale  will  be  developed. 

D.6.1  Development  of  Differential  Coefficients 

We  begin  by  casting  the  spatial  variation  of  a  function  in  the  operator  form 

6  =  6X  4-  +  50  ~  +  6H  ~  (0.25) 

0  A  0<P  Oil 

where  6X  and  6<f>  are  variations  in  longitude  and  latitude,  respectively,  and  are 
expressed  in  radians  while  5H  is  the  variation  in  height  expressed  in  units  of  the 
earth’s  radius. 

From  (D.3),  (D.5)  and  (D.6)  it  follows  that 

6u  =  —  6|,  6v  =  —  6t?,  6(h)  =  -  6X  (D.26) 

and  from  (C.25)  that 

jc_as_  =  0.  jc=as_=0  (d.27> 

0H  ~  0H  ’  0X 
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jC  =  9  sin:f>)  Vi  sin  cos  0  (D  28) 

3$  dip 


=  (l  -e:>  ~  -  SCJcJ  sin  0  cos  0 


(1X2  >) 


Consequently.  by  applying  (1X25)  to  £  ot  (1X1 )  wx  l,aVL‘ 


6£  -  A|  6X  +  A2  60  +  Aj6U 


(1X30) 


in  which 


A,  -  =  -  p  cos  0'  cos  @ 

dX 

A,  -  =  -  (SC5  +  H)  sin  0  sin  (h) 

3  30 

A3  =  =  cos  ^  sin  (h) 

oil 


(D.31 ) 


and,  similarly ,  for  v  of  (1X1) 


5t?  «  B.6X  +  Bj6<?  +  B3  6H 


(1)32) 


in  which 


B,  =  ill  =  _  {  sin  d 

'  ax 

B-  =  — 1  =  (SC:  +  H)  (cos  p  cos  d  +  sin  0  sin  d  cos  (h,  1 


B,  =  =  sin  Pcos  d  -  cos  0  sin  d  cos  ih; 

J  i  LI 


(1)33) 


However,  since  many  calculations  do  not  require  the  accuracy  of  (D.31)  and 
(D.33),  the  following  approximations  are  sometimes  used, 


A:  ~  ~  C5  p  sin  0'  sin  (h; 
A.,  *  p  cos  0'  sin  (h)  =  £ 


(0.34) 
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B,  ^  CS  p  cos  0'  cos  d  +  C2  p  sin  0'  sin  d  cos  (h) 


'  p  sin  0'  cos  d  -  p  cos  0'  sin  d  cos  (h)  =  rj 


(D.35) 


In  the  balance  of  this  sub-section,  those  differential  coefficients  associated 
with  6f.  6L;  6i  617,  6f,  6(fj2)  and  6n  will  be  developed  for  completeness. 

We  begin  by  applying  (D.25)  to  J  of  (D.l)  thereby  yielding 


6(  =  C,  5X  +  C2  60  +  C3  6H 


(D.36) 


in  which 


=  £  cos  d 


C2  -  =  (SC2  +  H)  (cos  0  sin  d  -  sin  0  cos  d  cos  (h)  ) 

00 

C3  =  -^  =  sin  0  sin  d  +  cos  0  cos  d  ccs  (h) 


(D  37) 


which,  like  the  approximations  of  (0.34)  and  (D.35),  can  be  approximated  as 
C2  =»  C  S  o  cos  0'  sin  d  -  C2  p  sin  0'  cos  d  cos  (h) 


C.  *=  p  sin  0'  sin  d  p  cos  0'  cos  d  cos  (h)  =  f 


(D.38) 


Further,  by  employing  (C.26)  or  (C.27)  which  define  Li  and  Lj ,  respectively,  we 
note  immediately  t!  at 


6L  =  -  tan  /  6f  =  -  tan  /  (C,  6X  +  C5  60  +  C3  6H] 


(D.39)* 


Next,  we  apply  (D.25)  to  £  given  by 


£  =  p  (C  +  H)  cos  0  cos  (h) 


*  is  usually  ignored  for  the  simple  reason  that  (  variations  can  enter  local  circumstances 
only  through  the  shadow  radii  and  thus  through  XL  of  (D.39).  Hence,  even  thoueft  0(£)  =  0(rj) 
-0(fi ,  tan  f  *=  0.C05  and  thus  6  L  is  sma>l  by  comparison. 
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and  find  that 


6£  =  A,'  6a  +  a;.  6tf>  +  A'j  6H 

in  which  (recalling  that  6  (h\  =  6X) 

A|  -  =  n  p  cos  <?  sin  (h  . 

3a 

A*2  =  =  -pi  (SC2  +H)sin  $  cos 

d(p 

A  3  =  =  ji  cos  <p  cos  (h 

Similarly,  (D.25)  applied  to  t?  of  (D.4)  yields 
5rj  =  Bj  6\  +  B2  6 <t>  +  B3  6H 

in  which 

B,  =  —  =  /isindA,  -  d  C, 

QK 

B2  =  =  {i  sin  d  Aj  -  d  C2 

o<p 

Bj  =  =  /i  sin  d  Aj  -d  C, 


0.40) 


(D.41 ) 


(D.42) 


(D.43) 


(D.44) 


(0.4.5) 


Hence,  we  can  write  immediately  (noting  that  6u  =  —  d £and  that  Sv  =  -  617) 

that 

6(nJ )  =  N,  6X  +  N2  60  +  N35H  (D.46) 

in  which 

N,  =  -  2(u  A',  +  v  B', ) 

Nj  =  -  2(u  A 2  +  v  B'2) 

N3  =  -2(uAi +vB'3)  (D.47) 

and  that 

5(n)  =  .  (D.48) 

2n 


D.6.2  Development  of  Differential  C^r  ctionsto  Loot1  Circumstances 

In  subsequent  calculations  of  the  differential  corrections  to  local  circum¬ 
stances,  8L,  fyi  =  -  8£,  5v  =  -  &rj,  6(n)and  6(n?)  will  be  neglected. 

Because  as  stated  previously,  the  penumbral  contact  tunes  are  not  required 
with  great  accuracy,  the  correction  will  be  based  on  (D.20)  and  thus 

6t  =  5  ( L*  -  m‘  ^  =  -  - 

\  2D  /  2D 

or,  since  near  the  contact  time  (L*  -  mJ)  =  0, 

6,  , _  AimD  =  u6*  +  v5t?  (D.50) 

2D  D 

Hence,  the  correction  to  the  penumbral  contact  time  can  be  written  as 

6t  =  PoX  +  q  6$  +  r  oil  (D.511 

in  which 

P  = 

q  = 


(uA,  +  vBj  )/D 
(nA2  +  vB2)/D 

(uA3  +  vB3)/D  (0.52) 


iLV  mI)5D  (D.49) 

2DJ 
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Similarly,  the  correction  to  the  time  of  greatest  phase  from  (D.14)  becomes 


»•„  '«(--£)  =  (P.S3) 

or 

6 1  =  p  6X  +  q  b<p  +  r  6H  (D.54) 

in  which 

Pm  =  (iiAi  +  vB,  )/n2 
qn,  =  (u  Aj  +  vB:)/n3 

rm  =  (u  A3  +vB3)/nJ  (D.55) 

The  correction  to  the  umbral  contact  times  will  be  calculated,  under  the 
assumption  that  (D.54)  has  been  calculated,  by  correcting  the  semi-duration  Su  of 
(D.20)  for  changes  in  locale.  This  is  accomplished  by  noting,  from  (D.19),  that 


and  that 


5K  =ps6X  +  qs54i  +  r$6H 

in  which 

Ps  =  (uB,  —  v  A,  )/n 
qs  =  (u  Bj  -  v  Aj)/n 
rs  =  (u  B3  -  v  Aj)/n 


Hence,  we  recalculate  the  semi-duration  Su  from  the  formula 


where 


cos  ^ 


sin  ru 


K  +  6K 
L2 


(D.57) 


(D.58) 


(D.59) 


(D.60) 
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Finally,  we  note  that  a  change  in  AT  from  AT  to  AT  +  6T  can  be  readily 
incorporated  within  the  existing  framework  by  noting  its  complete  equivalence  to 
a  fictitious  eliangc  in  longitude  (b<t>  =  6  H  =  0)  c.  the  form 

6\  =  1.002738  6T  (D.61) 

and  thus  can  be  accommodated  via  p,  pm  and  ps  (the  differential  coefficients  of 
5X)  in  computing  5t,  6tm  and  s'u  from  (D.50),  (D.54)  and  (D.57),  respectively. 
The  corrections  to  these  previous  times  (when  expressed  in  U.T.)  will  be  6t  -  6T 
and  6t  -  6T. 

m 

D.7  ECLIPSE  MAGNITUDE 

The  magnitude  of  an  eclipse  is  defined  as  the  fraction  of  the  solar  diameter 
covered  by  the  moon  at  the  time  of  greatest  phase.  Consequently,  referring  to 
Figure  D-l  —  a  redrawn  and  simplified  version  of  Figure  C-2  —  we  see  that  for  an 
observer  O^,  located  in  the  penumbral  shadow  PP'  that  the  magnitude  M,  of  the 
partial  eclipse  is  given  by 


which,  employing  the  similarity  of  triangles  PMOp  and  S'MA;  and  PMU  and  S  MS 
leads  to 


POP  _  PX-XOp 
PU  PX  -  XU 


But  PX  =  Li ,  the  penumbral  radius,  XOp=  m  and  XU  =  —  U2  during  totality  (as 
pictured),  hence 


M, 


L,  —  m 

L,  + \ 


(D.62) 


which  is  the  desired  result  valid  for  any  time  during  the  eclipse. 

Similarly  for  an  observer  located  within  the  umbra  at  Ou  the  magnitude 
is  given  by 


98 


r 
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But  B'B  =  SB  +  SB'  and  from  the  similar  triangles  UMOu  and  SMB;  and  PMOu  and 
S'MS  we  have  SB/S'S  =  UOu/PU  =  UOu/P'U,from  the  similar  triangles  SM'B'  and 
P'M'CV  and  S'M'S  and  U'M'P'  we  have  S'B/SS'  =  P'Ou/P'U',  hence 

SB  +  SB'  _  UOu  +  POu  _  P'U  _  P'X  +  XU 
SS  P'U'  P'U’  PX-XU'-  * 

or 

Mj  =  -  (D.63) 

L,  +  Li 


Alternatively,  we  can  write  that 


B'B 

S'S 


M’M  / 
— /Xy 


sin  St  Sw 

-  - 

sin  Sq  S0 


where  Sv  and  S0  are  the  semi-diameters  of  the  moon  and  sun,  respectively.  Hence 


^  *  M,  (D.64) 

S© 


0.8  DEGREE  OF  OBSCURATION 


The  degree  of  obscuration  at  a  particular  location  is  defined  as  that  fraction 
of  the  solar  disc  obscured  by  the  moon  during  the  partial  phase  and  is  —  in  the 
notation  of  Figure  D-2  and  based  upon  the  treatment  in  [2]  -  given  by 


S'= 


Ua  +aB) 


trr, 


© 


(D.6S) 


To  develop  S'  in  a  form  suitable  for  calculation  we  note  First  that 


r© 

>B  *2/ 

r0  cos  B 


^rO  -  x‘ 

J  dy  dx 

o 


or 
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which  becomes 


h  =  r© 


r©  cosB  v'tg  -  Xq  cos*B  +  rasin'"1  (cos  B) 


and,  since 


A  =  ~  -A,  , 


aB  =  r©  B  -  r0  sin  B  cos  B. 
Similarly,  we  can  write  that 

aA  =  r^A  -  r£  sin  A  cos  A 
and  hence  that  (D.55)  becomes 

S’  = 


^  A  +  B  -  (-7^“)  ^  A  cos  A  -  sin  B  00s  B 


r© 

Second,  we  note  from  (D.i  3)  that 

\ 

r© 


fD.66) 


(D.67) 


(D.68) 


S* 

V 


M2 


Hence,  if  we  set  the  solar  semi-diameter*  r@  =  BC  =  1 ,  we  can  write  -  referring 
once  again  to  Figure  D-2  - 

~~~  (Li  —  L5  ) 


r*  =  AC  =  s  =  — — 


(L,  +L,)  •  (D.69) 

From  (D.I  1)  we  note,  since  the  solar  remi-diameter  is  the  unit  of  distance,  that 


FD  =  2M,  =  - 


(L,  -  rti> 
(L,  +  L2) 


(D.70) 


and  thus  that  BA  1  BF  +  AF  »-  BD  —  FD  +  AF  -  1  —  2Mj  +  s,  or  that 


BA  = 


2m 


L,  +L; 


(D.71) 


1  It  should  be  noted  that  the  semi -diameter  of  the  sun  can  either  be  that  of  the  "usual”  optical 
sun  or  that  of  a  specific  radio  sun,  the  particular  semi-diameter  chosen  will  manifest  itself 
through  the  value  of  L2 


Next,  application  of  the  law  of  cosines  twice,  yields 


AC  1  =  BC  2  +  BA  2  -  2BC  BA  cos  B 


and 


BA2  =  AC 2  +  BC. 2  —  2AC  BC  cos  C 
or,  using (0.69),  D.70)  and  (0.71) 

cosB  =  (L,L2  +  m:)/'m  (L,  +  Lj)  (D.72) 

and 

oosC  =  (Lj  +  L  2-2inJ  )/(L,  +L})  (0.73) 

from  which 

A  ®  it  -  (B  +  C).  (0.74) 

Thus,  S'  of  (0.68)  becomes 

S' =  [s2  A  +  B  -  (s2  sin  A  cos  A  +  sin  B  cos  B)]/tt  (D.75) 

or,  since  the  third  term  is  twice  the  area  of  ABCA,  which  is  readily  shown  to  be 
'/i  s  sin  C, 


S'  =  (s2A  +  B  -  s  sin  C  ]n 


(D.76) 


where,  during  ihc  annular  phase.  S'  =  s3  and  during  the  total  phase  S'  =  1 


D.9  CENTROID  OF  UN0B9CURED  SOLAR  DISC 


From  symmetry  it  is  seen  that  the  position  of  the  unobscured  solar  disc 
centroid  lies  on  an  extension  of  the  line  joining  the  solar  and  lunar  centers  which, 
as  in  the  previous  section,  is  taken  as  the  x-axis  of  a  system  having  its  origin  at  the 
solar  center.  Hence,  we  may  write  -  referring  to  Figures  D-2  and  D-3  -  that  the 
centroid  of  area  B  (aB )  is  given  by 

vr^VrT3 

aB  ‘  x  b  =  2  f  ^  f  dy  x  dx 
r@cos  B  0 
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FIGURE  D— 3  CENTROID  GEOMETRY 
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or,  it  UK-  results  of  tiic  previous  section  are  employed 


3*  t  A  sin  A  cos  A) 
*  ( I  -  S') 


{  2nt  \ 

I  l.TlT  / 

'  / 


(D.83) 


Hence,  by  specifying  its  position  angle  =  l’  -rr  and  radial  distance  relative  to 
the  solat  center  xc  ,  the  centroid  may  be  located  easily  on  the  solar  disc. 

D  IO  TOPOCENTRIC  PARAMETERS 

The  topocentric  right  ascension  a ^  declination  6^  and  aistance  r^  of  the 
moon  can  be  calculated  from  the  following  exact  relationships 

r'k  cos  b ^  cos  =  rk  cos  6k  cos  <\  -  p  cos  0’  cos  r^s 
r^  cos  6^  in  cos  sin  ow  -  p  cos  0'  sin 

r^,  sin  6^  =  r^  sin  5 -  p  sin  0'  (D.84) 

in  which  crfc  and  are  the  geocentric  lunai  right  ascension  and  declination  of 
(A. 46),  respectively;  rv  is  the  geocentric  distance  given  by 

=  a^/sin  rrv  (D.85) 

where  the  lunar  parallax  rrc  <omcs  from  (A. 45),  and  rft  is  the  local  sidereal  time 
(L.S.T.)  at  r  hours  E.T.  given  by  (B.4). 

Next,  the  azimuth  Az^  (measured  east  from  north)  and  the  elevation  E£t  of 
the  moon  can  be  calculated  from 

cos  cos  Azw  =  +  sin  6^  cos  0  -  cos  5^  sin  0  cos  a'*) 

cos  sin  Az«,  -  -  cos  <5j.  sin  (r^  -  a^) 

sir.  =  +  sin  6^  sin  0  +  cos  cos  0  cos  (rt  -  a^)  (D.86) 

in  which  (rts  -  a^)  is  recognized  as  the  hoi,r  angle  of  the  moon. 
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And  finally,  the  topocentric  lunai  see, 
the  expression 


sir. 


S'fc  = 


sin  Sfc 


•erS’^,  can  be  calculated  from 


(D.87) 


where  the  sino  of  ti  e  geocentric  lunar  semi-diameter  can  be  rotten  directly 
from  (A.48)  we  note  in  passing  that  Sv  -  Sc  is  known  as  the  augmentation. 

In  the  case  of  the  sun,  an  approximate  form  of  < 0.84 )  suffices;  thus,  the 
topocentric  right  ascension  a©  and  the  declination  5©  are  given  by 

“©  =  a®  ~  ^Q© 

l  3  -  5©  -  A6©  (D.88) 


in  which 


Aa©  =  ttq  [p  cos  0'  sec  6©  sin  (TkS  -  a®) 

A60-  jr©  [p  sin  6'  cos  5©- p  cos  o'  sin  6©cos  (res  -  a©)]  (D.89) 

and  where  rr0  is  given  by  (A. 22).  As  in  CD. 86),  the  azimuth  Az0  and  KB©  are 
given  by 

cos  EB0  cos  Az0  =  +  sin  6©  cos  <5-  cos  6©  sin  0  cos  (r^  -  a©) 

cos  E2©  sin  Az  =  -  cos  6©  sin  (r^-  a'©) 

sin  EB©  =  +  sin  5©  sin  0  +  cos  6'©  cos  0  cos  (r^  -  aj©) 

(D.90) 

t\here  it  should  be  noted,  for  locations  near  the  center  line  and  at  or  around 
eclipse  maximum,  that  the  useful  approximation 

sin  E£0  —  £  CD. 9 1 ) 

follows  from  an  examination  of  (C.20). 

Since  the  solar  semi-diameter  is  almost  totally  insensitive  to  topocentric 
variations;  the  value  given  by  (A. 23)  can  be  used  throughout  the  eclipse. 
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APPENDIX  E 


SOLAR  COORDINATE  SYSTEMS  AND  RELATIONSHIPS 
E.O  INTRODUCTION 

In  this  appendix  various  coordinate  systems  applicable  to  solar  astronomy 
will  be  described  and  their  interrelationships  defined.  Particular  attention  will  be 
directed  to  the  development  of  both  precise  and  approximate  "pointing  instruc¬ 
tions”  for  observation  of  phenomena  both  on  and  above  the  solar  surface. 

E.1  PRELIMINARY  DISCUSSION  OF  SOLAR  COORDINATE  SYSTEMS 

Points  are  located  on  or  above  the  solar  surface  by  specifying  both  their 
heliographic  longitude  L  and  heliographic  latitude  B,  once  the  prime  meridian  of 
the  rotating"'  sun  (the  x-axis  of  the  heliographic  system)  and  the  rotation  axis  of 
the  sun  (the  z-axis  of  the  heliographic  system)  have  been  specified  for  the  date 
and  time  of  interest. 

This  specification  uses  as  reference  an  ecliptic  plane  having  no  secular 
motion  and  furthermore,  the  latitude  of  the  sun  is  neglected.  Hence,  one  requires, 
first,  the  inclination  of  the  solar  equator  with  respect  to  the  ecliptic  1,  or, 
equivalently,  the  angle  between  the  ecliptic  north  pole  zc  and  the  north  pole  of 
the  sun  zH  (this  is  shown  in  Figure  E-l);  second,  the  longitude  of  the  ascending 
node  of  the  solar  equator  on  the  ecliptic  U;  and  finally,  because  of  ihe  gaseous 
nature  of  the  sun  and,  hence,  the  lack  of  a  recognizable  reference  point  on  the 
solar  equator  with  which  to  measure  longitudes,  the  location  of  the  prime 
meridian  by  specification  of  the  heuographic  longitude  of  the  ascending  node  as  a 
given  function  of  time,  M.  These  ere  given  by 

1  =  7°15' 

o  =  7?°40'  +  50!*25  (t  -  !  850.0) 

M  -  180°  =  112?766  -  14°  184397 16  (J.D.  -  243  0000.5)  (E.l) 

in  which  t  represents  the  number  of  years  since  1850  and  J.D.  represents  the 
Julian  Day  number  of  the  date  and  time  of  observation. 

In  “actual”  appearance,  of  course,  the  sun  resembles  a  disc  onto  which 
points  both  on  and  above  the  solar  surface  are  apparently  projected.  The  normal 
to  the  disc  (x'  x'  in  Figure  E-l)  lies  along  the  earth-sun  line  at  the  celestial 
longitude  7T  where  is  given  by  (A.!  7)  o;  (AJ8).  A  projected  point  is 
located  on  the  disc  by  specifying  two  polar  coordinates:  first,  its  radial  distance 


*  In  what  follows  we  neglect  the  differential  {heliographic  latitude  dependent)  rotation  of  the 
sun.  Longitude  is  measured  westward  {in  the  direction  of  rotation)  from  the  prime  meridian 
and  ranges  from  3°  to  360°;  latitude  is  measured  from  the  solar  equator,  positive  towards  the 
north. 
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from  the  disc  center,  equivalent  to  the  angle  p  subtended  at  the  center  of  the  sun 
by  this  point  and  the  earth’s  center;  and  second,  its  position  angle  0  measured 
eastward  from  the  north  point  (the  z"  axis  of  Figure  li-1)  which  is  the  projection 
of  the  earth’s  rotation  axis  zL  on  the  disc*.  The  relationship  between  this 
“observational”  disc  coordinate  system  and  tire  heliographic  systems  can  be  fixed 
by  specification  of  the  position  angle  P  of  the  projection  of  the  solar  rotation  axis 
z(|  on  to  the  disc 

E.2  PRELIMINARY  CONSIDERATIONS 

In  the  remainder  of  this  appendix  we  assume  that  a  positive  rotation  with 
respect  to  an  axis  proceeds  >■  the  right  handed  sense,  and  furthermore,  in  the 
following  three  relations,  we  assume  that  the  rotation  operator  operates  upon  the 
coordinate  basis  x,  y,  z  to  pioduce  the  transformed  x’,  y\  /.'  basis. 


By  positive  rotations  with  respect  to  x,  y  and  z  axes  we  mean 


respectively. 


*  As  discussed  in  Appendix  D,  the  north  point  is  observalionally  defined  as  that  point  on  the 
solar  limb  which  is  the  tangent  point  of  the  parallel  of  declination  defined  by  the  apparent 
(diurnal!  motion  of  the  sun. 


Ill 


E.3  THE  HELIOGRAPHIC/ECLIPTIC  COORDINATE  TRANSFORMATION 


The  heliographi*;  system  (x„  ,  y„  ,  zH )  is  transformed  into  the  ecliptic  system 
(xt .  y, .  zv.  I  by  t Ih-  following  sequence  of  rotations:  first,  a  rotation  around  z()  by 
M;  second,  a  rotation  around  the  new  x-axis  by  -  I;  and  finally,  a  rotation  around 
the  new  z-.,xis  by  ft.  This  leads  to 


in  which 

Rc,,h  =  cos  ft  cos  M  +  sin  ft  sin  M  cos  1 
R1,]11  =  -  [sin  IT  cos  M  cos  1  —  cos  ft  sin  M  ] 

R*  j 1  =  sin  ft  sin  1 

R  j  ,H  =  sin  ft  cos  M  -  cos  ft  sin  M  cos  l 

R22 1  =  cos  ft  cos  M  cos  I  +  sin  ft  sin  M 

Ria11  =  -  cos  ft  sin  l 

Rji"  ~  -  sin  I  sin  M 

Rjj"  =  sin  1  cos  M 

Rj”  "  cos  I  (F..6) 

E.4  THE  ‘DISC '/HELIOGRAPH  1C  COORDINATE  TRANSFORMATION 
AND  RELATIONS 

The  ‘disc’/heliographic  transformation  can  be  developed  by  either  of  two 
routes,  namely:  that  the  ecliptic  system  (xt,  yt,  zc)  is  transformed  into  the  ‘disc’ 
system  (x‘,  y‘,  z‘)  by  a  rotation  around  zc  through  (Xq  +  7r)  and,  in  turn,  that 
the  ecliptic  system  is  transformed  by  use  of  (F..5)  or,  equivalently,  that  the  same 
development  as  that  leading  to  (E. 5)  is  used  but  with  -ft  replaced  by  -ft  +(Xq-ht  ) 
in  the  third  rotation.  In  either  case, 
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(E.7) 


with 


Rn  ~  _  cos  (XG  -  £2)  cos  M  +  sin  (Xq  -  £2)  sin  M  cos  I 
'H 

Rn  ~  —  [sin  (Xq  -  12)  cos  M  cos  I  +  cos  (Xq  -  £2)  sin  M  ] 

'H 

Ru  =  sin  (X@  -  £2)  sin  I 

f  h  __ 

R^  =  sin  (Xq-  £2)  cos  M  +  cos  (Xq  -  SI)  sin  M  cos  I 

#H 

R22  =  —  cos  (Xq  -  £2)  cos  M  cos  I  +  sin  (Xq  -  £2)  sin  M 

R»  =  cos  (Xq  -  £2)  sin  I 

Rj^  =  —  sin  I  sin  M 
Rjj  =  sin  I  cos  M 

R33  -  cos  I  (E.8) 

Next,  we  note  that  the  x'-axis  of  the  ‘disc’  system,  i.e.,  the  ‘disc’  center, 
defined  by  the  subterrestrial  point  located  at  heliographic  longitude  L0  and  lati¬ 
tude  B0 ,  can  be  specified  by  transposing  (E  7)  in  the  form 


or 


/  -  .  \ 
cos  a0  cos  Lq 

cos  B0  sin  Lo 

sin  B0 


=  (R  l1 


. 

\ 


/~'H  \ 
'  Ri:  ' 

R.2 


p'H 
1  3 


(E  9) 


cos  B0  cos  Lo  =  -  cos  (Xq  -  £2)  cos  M  +  sin  (Xq  -  £2)  sin  M  cos  I 
cos  B0  sin  Lo  -  -  [sin  (Xq  -  £2)  cos  M  cos  I  +  cos  (Xq  -  SI)  sin  M 1 

sin  B0  "  sin  (Xq  -  $2)  sin  I  (E  10) 
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These,  in  turn,  can  be  transformed,  by  means  of  some  stiaight forward  manipu¬ 
lation  into  the  more  commonly  employed  form  for  calculating  L0  and  B0, 

sin  Bfi  =  sin  (X0  -  LI)  sin  1 
cos  B0  cos  (l~o  —  N1)  =  ~  °°s  (\j)  -  ) 

cos  B0  sin  (Lo  -  M)  =  -  sin  <X0  -  LI)  cos  1.  (E.l  l) 


Finally  we  have  that  the  angle  Y  between  the  projection  onto  the  ‘disc’  of 
the  sun's  rotation  axis  z„  and  the  ecliptic  pole  can  be  gotten  in  two  steps:  first 
acquire  the  disc  coordinates  of  zlt  through 


(E.l  2) 


x'  =  sin  (Xq  -  fi)  sin  1 
v'  =  cos  (X0  -  iT)  sin  1 

z'  =  cos  I  (E.l 3) 


and  second,  note  that 

tan  Y  = - — 

2 


-  -  cos  (X0  -£2  )  tan  1. 


(E.14) 


E.5  THE  'DISCVEQ'JATORIAL  COORDINATE  TRANSFORMATION 
AND  RELATIONS 

The  ‘disc’/equatorial  transformation  is  developed  by  noting  that  the  equato¬ 
rial  system  is  transformed  into  the  ecliptic  system  by  a  rotation  around  xt  through 
the  obliquity  t  and  that,  in  turn,  the  ecliptic  system  is  transformed  into  the  disc 
system  by  a  rotation  around  zt  by  the  angle  (X0  +  it) .  Hence, 

A  / x  l- ' 

V  1=  IR'l  I  >r  1 

z  /  \z,  /  (E.l 5) 


/ 

l 
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with 


R.i 

- 

-  COS  Xq 

R'n 

- 

-  sin  X©  cos  c 

R’i'j 

= 

-  sin  Xq  sin  c 

K 

= 

sin  Xq 

Rjj 

= 

—  COS  Xq  COS  £ 

R» 

= 

-  cos  Xq  sin  e 

4 

= 

0 

R« 

= 

-  sin  c 

R33 

= 

COS  € 

i 


Furthermore,  we  also  note  that  the  angle  X  between  the  north  point  defining 
the  projection  of  the  earth’s  rotation  axis  zt  on  the  ‘disc’  and  the  ecliptic  pole 


can,  'ike  the  angle  \  be  gotten  in  two 
ze  through 


« 


or 


eps:  first,  acquire  the  disc  coordinates  of 
R'i'i  \ 

R&  I 

R33/  (E17) 


x’  =  -  sin  \q  sin  e 
y'  =  -  cos  Xq  sin  e 

z'  =  COS  £ 
and  second,  note  that 

tan  X  =  ~  -  -  cos  Xq  lane. 


(E.18) 


(F.19) 


j 
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Finally.  wc  see  tlut  the  position  angle  P  of  the  solar  rotation  axis  is  given  by 


P  -  X  +  Y 


(11-20) 


where  Y  was  defined  in  (1.14). 


E.6  THE  'DISCVDISC"  COORDINATE  TRANSFORMATION 


The  ‘discV'disc”  transformation  is  developed  by  rotating  the  ecliptic  pole  (the 
z'-axis)  around  the  x'-axis  by  the  angle  -  X  into  the  north  point  (the  /'-axis), or  that 


I  0  0 

1  cos  X  -  sin  X 
0  sin  X  cos  X 


(F..21 ) 


E.7  THE  HELlOGi\APHtC/"DISC  '  COORDINATE  TRANSFORMATION 
AND  RELATIONS 


This  "disc'’  system  (x',-  y"  z")  is  transformed  into  the  heliograpliic  system 
(x„  ,  yH  .  zH  )  l  y  the  following  rotations  first.a  rotation  around  x"  by  P;  second,  a 
rotation  around  the  new  y-axis  by  B0;  and  finally,  a  rotation  around  the  zH  axis 
by  -  Lo .  This  leads  to 


(E.22) 


in  which 

R,V  -  c\)s  B0  cos  1  o 

RiV  =  -  (cos  P  sin  Lq  -  sin  P  cos  Lo  sin  B0  1 

R’,V  =  -  l^in  P  sin  Lo  +  cos  P  cos  Lo  sin  B0  ) 

R-,1,"  -  cos  Bq  sin  L0 

R'jV  =  cos  P  cos  L0  +  sin  P  sin  L0  sin  B0 

R'A"  -  sin  P  cos  L0  -•  cos  P  sir.  sin  B0 

KVi  =  sin  B0 

R'jV’  =  -  sin  P  cos  B0 

rV3  -  cos  P  cos  B0 


1  l(. 


(H  23) 


a fd  longitude  L. 


ily,  a  point  «.-i  the  surface  of  the  sun  (of  radius 
,  tiKK'by  having  tlu  heliographic  coordinates 


R©  )  at  latitude  H 


-v  | ■  ■*  R  ©  cos  If  cos  L 
'  ll  =  R©  cos  B  sin  l 
z„  =  R©  sin  1 i. 


n  24  > 


Will  be  transformed  into  a  point  in  the  “disc"  system  with  heliocentric  angle  p 
and  position  angle  0  thereby  having  the  “disc”  coordinates 

x  ~  R©  cos  p 

-  -  R©  sin  p  sin  0 

z"  ~  R©sinpcosd  (|; ->5) 

5f  thes<  arc  substituted  in  0.22),  and  appropriate  operations  and  manipulation' 
performed,  we  have  that 


sin  B  -  sin  B0  cos  p  +  cos  B0  sin  p  cos  (P  -  0)  d  26) 

cos  B  sin  (L  -  ].„)=■  sin  p  sin  ( ( F  -  0)  d;  -’?) 

which  permits  the  immediate  calculation  of  the  heliographic  latitude  and  longi 
tude  of  a  point  once  its  disc  coordinates  p  (whose  determination  will  be  discussed 
in  Section  £.10)  and  0  have  been  specified.  Equivalent!  •,  we  have 

cos  p  =  cos  B0  cos  B  cos  <  L  -  L0  )  +  sin  B0  sin  B  (f  .2«) 

sin  p  sin  (p  ~  0)  ^  cos  B  sin  i  L  -  L„  )  d-  ->0( 

if  B  and  L  are  specirieel  and  p  and  0  are  equired 

E.8  RIGHT  ASCENSION  AND  DEC Li.vATlUd  OF  AN  OBSERVED  DOlNT 
ON  THE  "DI9C"  . ‘ 


This  development  is  initiated  by  employing  the  transpose  of  ( I 


on  the  ‘Hjs.- 


.  ►!,„  . .  ._  ...  .  . 

mv  i^/mi  ill  lllc  lOI  111 


1 )  operating 


x'  =  R©  cos  p 

v'  “  R©  sin  p  sm  (G  - 

7  ~  ~  R©  sin  p  cos  (0  -  X)  d  fOi 


'17 


to  yield 


xi  ~  R<$>  I  _  cos P  cos  Xq  -  si np  sin  X^sin  (0  -•  X)] 

y,  =  R(j,  I  -  cos  p  sin  Xq  cos  e  +  sin  p  cos  Xq  cos  e  sin  (0  -  X) 

-  sin  p  sin  e  cos  (0  -  X)) 

z,  =  Rq  1  -  cos  p  sin  Xq  sin  e  +  sin  p  cos  Xq  sin  e  sin  (0  -  X) 

♦  sin  p  cos  e  cos  (0  -  X)] .  (E.31) 


Next,  because  the  geocentric,  equatorial  x,  y,  z  components  of  the  solar  center 
(neglecting  the  solar  latitude)  are  given  by  (A.21)  as 

Xq  =  R  COS  6 q  COS  Oq  =  R  cos  Xq 
y  0  =  R  cos  6 0  sin  Oq  =  R  sin  Xq  cos  e 

Z0  =R  sin  6q  =  R  sin  Xq  sin  e  (E.32) 


then  the  geocentric,  equatorial  components  of  tne  disc  point  can  be  written  by 
x  =  X0  +  xF  =  R 


1  -  cos  P  )cos  X- 

R  / 


sin  P  1  sin  (0  -  X)  sin  Xq 


y  =  y®  +  yr  "  R 


/, _  fa  .V 

j  ^  (>  .’5 ,J  I  sin 

\  R  / 


Aq  cos  e 


+  f  sin  pjtsin  (0  -  X)  cos  e  cos  Xq  -  cos  (0  -  X)  sin  e ) 


z  =  Z0  +  zL  =  R 


1 - ^  cos  Pj  sin  Xq  sir;  t  + 


+  | — —  sin  p  |  (cos  (0  -  X)  cos  e  +  sin  (0  -  X)  sin  e  cos  Xq) 


(E.33) 
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Conseq'iently,  if  we  employ  the  reasonable  approximation  that  x2  +  y 1  ♦  r 
R5 ,  we  cart  write  immediately  that  the  nght  ascension  at  and  declination  b%  of 
the  point  are  given  by 

R  cos  Ss  cos  as  =  x 
R  cos  i .  sin  as  =  y 

R  sin  Ss  ~  z  (F..34) 

where  x,  y,  and  z  are  given  by  (E.33). 

But  (E.34)  and  (E.33)  can  be  transformed  into  a  more  familiar  form  by 
expanding  the  sin  (0  —  X)  and  cos  (0  -  X)  terms  and  identifying  sin  x  and  cos  x 
with  the  ecliptic  and  equatorial  coordinates  of  the  sun.  For  from  (E.  1 8)  and  (E.32) 
we  l lave  that 

(y')a  +  (z')3  =  (-  cos  Xq  sine  )2  +  cos2  e  =  cos2  (E.35) 


so  that 


sin  X  = 


y* 

VfTFTTi7? 


—  cos  X  0  sin  e  /cos  6 0 


(E.36) 


and 


cos  X  = 


cos  e./cos  U0 


(E.37) 


'^(y't3  -Hz')' 


which,  when  employed  -  together  with  (E.32)  -  in  (E.34),  yield  (after  consider¬ 
able  manipulation) 


cos 


6  cos  a  =  / 1  -  cos  p )  cos  50  cos  a0 
s  5  V  R  / 


R, 


sin  p  (sin  0  sin  a^+  cos  0  cos  a0  sin  60] 


/  Rra  \ 

cos  6S  sin  a  =  <  1 - ^  cos  p  j  cos  60  sin  o0 


+ 


sin 


sin  p  [sin  0  cos  a0-  cos  0  sin  a0  sin  6 q  ! 

(R  \  R 

1 - —  cos  p )  sin  6/7)  +  — —  sin  p  cos  0  cos  6,7 1  , 

R  /  u  R 


( 12.38) 
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K0 

— ^  s>n  p  sin  0  -  sin  (u^  -  a0)  cos  5^ 

R 


V-i 


sin  o  cos  0  -  cos  60  sin  6s  -  sin  60  cos  6s  cos  (ov  -  a0) 


(  R0  \ 

^  1  -  cosp|  =  sin  60  si 


sin  6s  +  cos  3t-  cos  5s  cos  (a,  -  a0).  (E.39) 

But  reference  to  Figure  H-2  shows  tlwt 


r< 


3 


R, 


« 


and 


sin  p,  -  sin  P  *=  — ~  sin  p 

u?  s  K 


R  -  R0  cos  P  R0 

cos  p  |  -  - — -  %  1 - —  cos  p 

«t  R 


(E.4C) 


(E.41) 


so  that  (E.39)  reduces  to  the  fainihar  set  of  relations 
sin  Pj  sin  0  -  sin  (os  -  a ^  cos  6$ 

sin  p,  cos  6  =  cos  6@  cos  6s  -  sin  60  sin  6s  cos  (a  a  0 ^ 

cos  p,  -  sin  60  sin  6s  +  cos60cos  5s  cos(as-OQ  )  (F.42) 

which  also  follow  directly  from  Figure  E-3  and  the  solution  of  the  astro*  mica! 
triangle  on  the  celestial  sphere  of  infinite  radius  (wherein  the  approxirru  .ons  of 
(E.40)  and  (E.41)  are  rendered  unnecessary).* 

If  the  separation  of  (a0,  50)  and  (a,. .  6s)  is  sufficiently  small  that  square  and 
higher  order  terms  can  be  neglected,  then 

p,  sin  0  (o*s  -  a0)  cos  6^ 

p,  cos#  *  sin  (fis  -  60)  '  &s  --  fi0  , 


*  See  [91  pp.  "15-26. 
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or,  to  withm  tins  level  ol  approximation 

(t.43) 
(t  44) 

E.9  DETERMINATION  OF  p  FROM  OBSERVATIONS 


Qy  -  a0  *  p  i  sin  0  see  60 
-  60  ==  p,  cos  0 . 


The  determination  of  p  is  based  on  the  geometry  of  Figure  F.-4,  a  more 
detailed  version  of  Figure  li-2,  which  also  provides  for  a  radio  sun  of  radius 
Rr  =  3fR  R0  with  XR  >  1  as  well  as  for  the  optical  sun  of  radius  R.  r  and  r0 
are  meant  to  symbolize  linear  observational  measures  in  some  convenient  scale  on 
say, a  photograph. 


Because  of  its  small  magnitude  16')  we  can  write  that  the  solar  semi¬ 
diameter  of  A. 2 3  is  given  by 


sin  S0  -  —  *  S © 


(E.45) 


and,  to  the  same  level  of  approximation  that 


Pj_ 

se 


r 


(E.46) 


Hence,  p,  can  be  determined  from  the  three  known  or  measured  quantities  S^,, 
r  and  r0 . 


But  we  see  that 

s*n  A  _  sin  (p  |  +  p) 

R  R 

or,  employing  (E.45),  that 

sin  ip,  +  p)  *  Pi /S0  (E.48) 

where  it  should  be  pointed  out  that  the  optical  radius  has  bee  .  tacitly  employed. 

If  on  the  other  hand  we  wish  to  calculate  the  p  appropriate  to  a  radio  sun  of 
radius  Rr  =  JfR  R0  from  the  same  p, ,  we  modify  (E.48)  such  that 

«n<p1+p,l  •  P./S**  CE  W 


sin  p, 

R/n 


£l 

R© 


(E.47) 
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